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PREFACE 


This  Is  one  of  a series  of  Rand  reports  documenting  a study  of 
advanced  employment  concepts  for  ground  force  operations,  sponsored 
hv  the  Tactical  Technology  Office  of  the  Defense  Advanced  Research 
Projects  Agency  i ARPA) . The  overall  objective  of  the  study  was  to 
develop  new  concepts  for  the  employment  of  ground  combat  systems  in- 
corporat ing  advanced  technology,  and  to  define  and  evaluate  weapon 
systems  for  implementing  these  concepts. 

The  first  report,  R-23t>5-ARFA,  A Method  for  Evaluating  Advanced 
Susteme  aid  Concept*  for  (Around  Combat,  by  E.  W.  Paxson  and  M.  G. 
Weiner,  describes  the  general  method  of  the  study  and  discusses  spe- 
cittc  stops  with  examples  from  the  evaluations  carried  out  to  date. 

The  method  incorporates  a three-dimensional  terrain  board,  a computer 
model,  and  a series  of  analytic  modules  programmed  for  a hand  calcula- 
tor to  assess  the  outcomes  of  various  combat  engagements. 

The  present  report  describes  in  detail  the  computer  model,  TIMER, 
used  in  the  evaluation  method  for  determining  the  et feels  ot  terrain 
on  visibility.  It  is  being  published  separately  for  the  convenience 
of  readers  who  will  find  TIMER  useful  in  applications  other  than  t lie 
one  described  here. 

A third  volume,  R-2377-ARPA,  Interactions  bW:.\vu  Tactics  , ptd 
Technology  in  Cnumd  W,irJ\iro,  by  M.  G.  Weiner,  E.  W.  Paxson,  and 
R.  A.  Wise,  presents  a series  of  examples  of  interactions  between 
tactics  and  technology  that  were  observed  in  the  combat  evaluations 
conducted  during  the  study.  These  examples  form  the  basis  for  some 


general  observations  and  speculations  about  tactical-technological 
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SUMMARY 


The  Terrain  Intervisibility  and  Movement  Evaluation  Routine  (TIMER) 
was  developed  as  part  of  a Rand  study  of  the  relationships  between  the 
tactical  employment  and  the  technological  capabilities  of  advanced  weap- 
on systems  for  ground  combat.  The  purpose  of  TIMER  was  to  determine 
the  effects  of  terrain  on  target  visibility  during  combat  engagements 
carried  out  on  a terrain  board  by  human  players.  Among  the  questions 
that  TIMER  was  designed  to  supply  data  on  were:  the  occurrence  and 
duration  of  target  visibility  (line  of  sight)  from  defense  positions, 
the  maximum  range  at  which  target  acquisition  might  occur,  and  the  po- 
tential rate  at  which  a defense  force  could  fire  on  ("service")  an 
attacking  force. 

The  TIMER  model  uses  the  digitalized  terrain  data  base  developed 
by  the  Defense  Mapping  Agency.  With  this  fairly  fine-grained  data 
base,  the  program  is  able  to  determine  whether  there  is  a line  of  sight 
between  any  two  points  on  the  terrain.  From  this  basic  computation, 
a variety  of  conditions  can  be  represented.  It  is  possible  to: 

o Input  one  or  more  potential  avenues  or  routes  of  attack  and 
calculate  the  visibility  of  targets  on  the  routes  from 
various  defense  positions. 

o Specify  different  target  velocities  and  calculate  the 
length  of  exposure  of  targets  on  the  routes. 

o Incorporate  the  reaction  time  of  different  defense  systems 
and  determine  the  number  of  opportunities  to  engage  targets. 

o Determine,  for  a given  defense  force  occupying  a number 

of  different  defense  positions,  their  potential  ability  to 
"service"  an  attacker  force  approaching  on  one  or  several 
avenues  of  advance. 

The  report  includes  the  TIMER  inputs,  program  listings,  and  out- 
puts for  these  conditions,  and  furnishes  a number  of  illustrations  of 
the  results  that  TIMER  can  produce. 
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I.  GENERAL  DESCRIPTION  OF  TIMER 


j 

The  Terrain  Intervisibility  and  Movement  Evaluation  Routine  (TIMER) 
is  a computer  model  that  was  developed  as  part  of  a study  of  advanced 
employment  concepts  sponsored  by  the  Tactical  Technology  Office  of  the 
Defense  Advanced  Research  Projects  Agency  (ARPA) . This  study  examined 
the  relationships  between  the  tactical  employment  and  the  technological 
capabilities  of  advanced  weapon  systems  for  ground  combat. 

A method  for  evaluating  those  systems  was  devised;  it  uses  a three- 
dimensional  terrain  board,  human  players,  and  computer  and  analytic 
programs  for  conducting  a minute-by-minute  evaluation  of  a combat  situa- 
tion involving  the  advanced  weapon  systems.  TIMER  was  developed  to 
produce  data  on  the  effects  of  terrain  on  target  visibility;  among  the 
questions  it  was  designed  to  answer  were:  the  occurrence  and  duration 
of  target  visibility  (line  of  sight)  from  defense  positions  in  specific 
terrain,  the  maximum  range  at  which  target  acquisition  might  occur,  and 
che  potential  rate  at  which  a defense  force  could  fire  on  ("service") 
an  attacking  force. 

The  current  version  of  TIMER  uses  European  terrain  in  the  border 
area  between  the  Federal  Republic  of  Germany  and  the  German  Democratic 

2 

Republic  in  the  U.S.  V Corps  sector.  An  area  of  approximately  2000  km 
is  represented.  The  terrain  data  base  for  the  area  consists  of  digi- 
talized data  obtained  from  the  Defense  Mapping  Agency,  and  includes  a 
data  point  at  every  12.5  meters  in  the  area.  For  each  data  point,  the 

elevation  to  the  nearest  meter  is  given,  as  well  as  whether  the  point 

* 

is  in  a forest  or  urban  area. 

With  this  fine-grained  data  base,  the  TIMER  program  permits  the 
user  to  specify  one  or  more  defense  position  locations,  down  to  the 
level  of  individual  vehicles  as  well  as  one  or  more  specific  targets. 

The  program  determines  whether  there  is  a line  of  sight  between  each 


The  Defense  Mapping  Agency  digitalized  data  base  for  the  Federal 
Republic  of  Germany  covers  a larger  area  than  that  used  in  the  TIMER 
program,  but  the  current  program  can  be  upgraded  to  a larger  area. 


r 
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defense  position  and  the  targets.  From  this  basic  computational  capa- 
bility, it  Is  possible  to: 

o Input  one  or  more  potential  avenues  or  routes  of  attack 
and  calculate  the  visibility  of  targets  on  these  routes 
from  a variety  of  different  defense  positions. 

o Specify  different  target  velocities  and  determine  the 
length  of  exposure  of  the  targets  on  the  routes. 

o Incorporate  the  reaction  times  of  different  defense 
systems  and  determine  the  number  of  opportunities  to 
engage  enemy  targets. 

o Determine,  for  a given  defense  force  occupying  a number 
of  different  defense  positions,  their  potential  ability 
to  "service"  an  attacker  force  approaching  on  one  or 
several  avenues  of  advance. 

These  and  other  conditions  and  situations  can  be  incorporated  in 
the  model  to  provide  quantitative  data  on  the  effects  ot  terrain, 
vegetation,  and  urban  build-up  of  a geographic  area  on  the  potential 


of  combat  systems. 
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II.  INTERVISIBILITY  CALCULATIONS 

The  fine-grained  TIMER  data  base  described  in  the  Introduction  and 
a basic  point-to-point,  line-of-sight  calculation  have  been  incorporated 
in  the  TIMER  computer  program.  It  determines  intervisibility  between 
defense  positions  and  attacker  vehicles  moving  along  specified  paths 
or  routes.  Appendix  A presents  a description  of  the  data  base  and  a 
program  listing  for  the  intervisibility  computations. 

The  steps  involved  in  using  the  program  include: 

o Selection  of  routes  or  paths , which  requires  detailed 
examination  of  topographic  maps.  Ancillary  data  not 
always  found  on  maps  may  also  be  needed,  such  as  the 
widths  and  load-bearing  capacity  of  roads  through  forested 
areas.  The  selection  of  preferred  routes  must  take  into 
account  not  only  trafficability,  cover,  and  concealment 
afforded  by  the  various  topographic  features,  but  also 
expected  tactics  employed  by  enemy  mechanized  units. 

Each  route  selected  is  defined  for  the  computer  as  a 
sequence  of  points  about  200  meters  apart,  and  the  route 
is  taken  to  be  linear  between  consecutive  points. 

o Selection  of  defensive  positions , which  also  requires 

detailed  examination  of  maps  to  take  topographic  details 
into  account  plus  such  factors  as  expected  attack  routes, 
number  of  defense  positions,  availability  of  withdrawal 
routes  for  the  defenders,  and  characteristics  of  the 
defender  and  attacker  weapons  systems — maximum  range, 

' ' 

reaction  time,  velocity,  and  the  like. 

Figure  1 illustrates  a set  of  attack  routes  and 
defense  positions  selected  in  a study  involving  two  ad- 
vanced systems,  a laser  beam-rider  missile  system  for 
direct  fire  and  an  indirect  fire  weapon  system  that  ac- 
quired targets  by  means  of  an  elevated  sensor.  The  units 
using  these  systems  occupied  defense  positions,  indicated 

I I 

_ _ I: 


• Location  of  defense  unit 
Major  attack  routes 
«■*  Secondary  attack  routes 

Fig.  1 — Initial  defense  force  deployment  and  enemy  attack  routes 
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by  dots  in  Fig.  1,  to  provide  initial  coverage  of  the 
likely  (but  unknown)  routes  of  enemy  advance.  In  one 
hypothetical  situation,  the  enemy  attacked  with  two  regi- 
ments forward,  and  each  regiment  fanned  out  to  place  one 
battalion  on  each  of  three  routes.  The  heavy  solid  lines 
in  the  figure  indicate  the  two  major  attack  routes,  pri- 
marily using  roads;  the  heavy  dashed  lines  indicate  four 
secondary  attack  routes,  two  flanking  each  major  attack 
route. 

The  coordinates  of  each  defense  location  are  input  to  the  computer 
program  along  with  the  following  parameters  that  further  characterize 
or  constrain  the  capability  of  the  individual  observers: 


Ranges  of  observers/sensors  (maximum  and  minimum,  km) 
Coverage  sector  (angles  measured  from  north) 

Observer  height  (meters) 

Target  height  (meters) 

Forest  and  urban  height  (meters) 


The  computer  program  determines  intervisibility  between  each  defense 
position  and  each  attack  route  for  points  along  the  route  (at  the  input 
height  above  the  route)  a fixed  increment  apart  (e.g.,  25  m,  an  input 
parameter)  taking  into  account  the  minimum  and  maximum  range,  coverage 
sector,  observer  and  target  heights,  and  forest  and  urban  heights. 

Table  1 illustrates  tne  TIMER  output  for  a single  position  on  a single 

attack  route.  The  table  presents  the  route  segment  intervisibility  in 

* 

a condensed  form. 

The  input  values  for  the  TIMER  run  are  listed  at  the  top  of  the 
table.  The  coordinates  of  the  observation  point  are  65.20  km  Easting 
and  121.70  km  Northing,  and  the  viewing  sector  is  between  30  and  150 
degrees  measured  from  north.  The  surface  elevation  at  the  observer 
location  is  obtained  from  the  data  base.  The  0 value  for  "observer 


The  output  from  TIMER  all j includes  a disk  file  data  set  contain- 
ing the  range  to  each  visible  segment  on  each  route. 
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Table  1 

VIS1BI.K  RODTK  STRKTCHKS 


Parameters 

Observer  location: 

X - 65.20  km  Y - 121.70  km 

Observer  symbol: 

K92 

Route  symbol: 

C4 

Viewing  sector: 

10  to  150  degrees  clockwise  from  north 

Minimum  range: 

0.75  km 

Observer  height: 

30  m 

Target  height: 

2 m 

Forest  he Ight : 

15  m 

Town  height: 

10  m 

t'alculated  Values 

Observer  location  surface  elevation:  278  m 

Observer  location  topography:  0 


Rout  e 

Maximum  Range 

Seg. 

Xcoord 

Ycoo  ril 

nist 

Hear 

0.75 

l . 50 

1.0 

101 

66.4‘> 

22.60 

1.6 

55.  1 

0 

0 

l 

10  1 

66 . 44 

22.59 

1.5 

54.3 

0 

0 

3 

125 

65 . 94 

22.41 

1.0 

45.3 

0 

l 

1 

l 15 

65.  72 

22.27 

0.8 

42.5 

0 

l 1 

1 1 

location  topography"  Indicates  that  the  observer  location  Is  nol  In 
a forest  or  town. 

The  bottom  part  of  the  table  shows  the  output.  It  lists  the  route 
segment,  coordinates,  distance,  and  bearing  from  the  observer  or  sensor 
location  at  which  Intorvislbil  Ity  starts  and  ends  tor  the  three  sensor 
ranges  given.  For  example,  at  the  maximum  sensor  range  ot  J.O  km, 
Intorvislbil  Ity  first  occurs  at  segment  1.01  (the  101st  25-m  segment 
from  the  start  o>  the  route)  and  lasts  to  segment  103,  a total  ol  three 
segments,  or  a stretch  of  75  m.  lnterv.lsibll.ity  does  not  occur  again 
until  segment  125  but  then  continues  until  segment  135,  a total  ot  ll 
segments . 

For  a maximum  sensor  range  of  1.50  km,  intorvislbil Ity  occurs  only 
for  segments  125  to  l J S , because  segments  101  to  103  are  beyond  the 
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1.5  kiu  range.  And  for  a maximum  sensor  ran  no  of  0.75  km,  no  intervisi- 
bllities  are  calculated,  because  the  minimum  range  specified  in  the 
Inputs  is  0.75  km. 

A similar  calculation  is  made  from  every  observer  or  sensor  loca- 
tion to  every  route  within  the  view  lug  sector  specified  in  the  input. 
However,  it  does  not  present,  in  a conventional  form,  an  Integrated 
picture  of  the  intervisibility  from  all  of  the  observation  points  on 
all  ol  the  routes.  This  is  obtained  by  combining  a TIMER  output  data 
set  containing  the  ranges  to  all  of  the  visible  segments  with  a modified 
TIMER  run  indicating  which  segments  are  not  visible  because  they  are  in 
a forest  or  town.  Figures  2 and  J illustrate  the  results  ol  this  pro- 
cedure. Figure  2 shows  the  stretches  of  the  six  attack  routes  in  Fig.  1 
that  are  in  the  "open"  and  the  stretches  that  are  in  lorests  and  towns. 
Only  about  half  the  total  length  of  the  routes  is  in  the  open,  although 
individual  routes  vary  from  22  percent  to  82  percent  open.  l’he  "open" 
areas  represent  /;/  visible  portions  of  the  routes,  but  not 

all  of  them  are  actually  visible  from  the  defense  positions,  because 
hills,  forests,  and  towns  may  interfere  with  the  line  ol  sight.  Fig- 
ure l shows  the  percentages  of  the  routes  that  are  actually  visible 
forward  from  one  or  more  of  the  defense  positions,  for  a situation  in 
which  tlie  sensors  of  the  indirect  fire  systems  are  elevated  to  a height 
of  JO  meters.  About  JO  percent  of  all  the  routes  combined  are  "covered" 
by  this  elevated  sensor,  and  about  50  percent  ol  the  potential Iv  visible 
stretches  are  covered. 

The  examples  presented  in  this  section  Indicate  how  t lie  TIMER 
model  can  be  used  to  determine  where  and  how  much  coverage  ol  potential 
attack  routes  can  In'  obtained  from  one  or  more  defense  posit  tons  with 
elevated  sensors.  Although  not  illustrated  here,  the  coverage  I rom  each 
defense  position  is  obtained,  as  well  as  the  amount  and  local  ion  ot 
overlapping  coverage  from  other  defense  positions.  Similar  data  can 
be  obtained  for  defense  positions  that  do  not  have  elevated  sensor 
plat  I orms. 


Percent 


8 8 8 5 8 
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111.  FIRING  OPPORTUNITY  CALCULATIONS 


As  described  in  Sec.  I l,  the  TIMER  program  calculates  the  locations 
and  lengths  of  all  the  visible  stretches  on  each  attack  route  for  each 
defense  position,  within  the  constraints  Imposed  upon  the  observer. 

As  Illustrated  in  Fig.  3,  the  number  and  lengths  of  the  visible  stretche 
can  varv  considerably  both  along  the  route  and  from  route  to  route. 

The  extent  of  the  Intervisibility  of  a particular  defense-observer/route 
combination  can  be  represented  by  such  metrics  as  the  number  and  average 
length  of  the  visible  stretches,  but  intervisibility  data  should  be 
related  directly  to  defense  system  characteristics.  To  this  end,  a 
metric  called  a "firing  opportunity"  was  developed  during  the  course 
of  the  study  and  introduced  into  the  TIMER  program.  The  firing  oppor- 
tunity takes  two  factors  into  account: 

o The  length  of  time  an  attack  vehicle  Is  exposed  to  fire 
on  a visible  stretch  depends  on  the  speed  of  the  vehicle. 

o The  ability  to  deliver  ordnance  on  the  attack  vehicle 
from  the  defense  positions  depends  on  how  long  it  takes 
the  defense  system  to  respond,  once  a target  is  in  view. 

This  system  "reaction  time"  is  defined  as  the  time  from 
first  line-of-sight  contact  on  a visible  stretch  until 
the  ordnance  arrives  at  target.  That  Is,  it  Includes 
sighting,  launching  and  time  of  flight  to  the  target. 

For  example.  If  a vehicle  is  moving  at  30  km/hr  (300  m/min)  and 
tin1  system  reaction  time  is  two  minutes,  the  system  requires  a visible 
stretch  of  1000  in  from  initial  line  of  sight  to  final  ordnance  delivery. 
Each  1000  m of  visible  stretch  thus  provides  one  firing  opportunity. 
Processing  the  computer  output  data  on  visible  stretches  makes  it  pos- 
sible to  determine  the  number  ot  firing  opportunities  for  different 
cond  it  Ions . 

Figure  4 illustrates  how  the  number  ol  firing  opportunities  can 
be  used  in  evaluating  the  tradeoff  between  system  character  1st Ics.  In 


Firing  opportunities 

Fig.  4 — Guided  mortar  firing  opportunities  for 
various  reaction  times 

the  figure,  tiring  opportunities  are  given  for  an  indirect  tire  system 
that  lias  an  elevated  sensor  platform  for  target  acquisit ion.  The  t ig- 
ure  shows  total  firing  opportunities  as  a function  ot  platform  height 
and  system  reaction  time.  As  expected,  firing  opportunities  increase 
as  Lite  height  of  the  sensor  increases  or  as  the  system  react  ion  time 
decreases.  The  results  provide  one  basis  for  determining  the  value  ot 
either  increasing  the  sensor  height  or  decreasing  the  system  reaction 
time,  and  also  the  tradeoff  between  the  two.  For  example,  it  tin'  nonti 
nal  platform  height  is  50  m and  the  nominal  reaction  time  is  1 JO  sec, 
doubling  the  platform  height  to  100  m doubles  the  number  ot  tiring 
opportunities  from  5 to  10,  whereas  halving  the  sytem  reaction  time 
from  1J0  see  to  ftO  sec  more  than  triples  the  number  ot  tiring  oppor- 
tunities, from  5 to  In. 
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IV.  SKTV  U:  INy  KiVrK  yAI.yjll._AT IONS 


The  "1'lrlnR  opportunity"  metric  do  lined  In  Sec.  Ill  Is  a measure 
of  the  potential  of  an  Individual  defense  position  against  an  individual 
attack  route.  To  measure  the  combined  capabilities  of  a set  of  defense 
positions  against  a set  of  attack  routes,  we  total  the  firing  opportuni- 
ties from  all  defense  positions  against  all  the  attack  routes.  However, 
certain  factors  could  diminish  the  number  of  firing  opportunities  on 
enemy  units  advancing  along  the  routes  in  a coordinated  attack.  For 
example,  some  of  the  routes  may  be  much  less  visible  than  others,  and 
the  enemy  may  choose  to  concentrate  his  forces  on  those  routes,  or  the 
timing  of  the  attack  may  saturate  the  defenders'  capabilities  to  use 
all  the  tiring  opportunities  presented. 

To  account  for  some  of  these  factors,  a computer  routine  was  added 
to  the  basic  TIMER  routine  as  part  of  a study  done  for  the  II. S.  V Corps 
and  the  results  presented  in  terms  of  a new  metric,  the  "servicing  rate" 
of  the  defending  units.  The  servicing  rate  may  be  defined  by  a simple 
example:  A defending  force  that  fires  upon  15  tanks  in  live  minutes 
wll'  he  said  to  have  a servicing  rate  of  three  tanks  per  minute.  Al- 
though the  servicing  rate  of  a given  force  is  situation-dependent  and 
the  computer  model  Is  a simulation  of  a particular  situation,  the  illus- 
tration presented  below  covers  a number  of  aspects. 

The  situation  is  one  in  which  enemy  tank  companies  consisting  of 
ten  tanks  each  are  attacking  along  each  of  ten  routes  with  an  approxi- 
mate spacing  of  50  m between  tanks.  All  companies  are  moving  at  the 
same  constant  speed  of  18  km/hr  or  5 m/sec.  The  starting  positions  on 
the  routes  are  inputs  and  may  be  arbitrary,  but  they  may  be  chosen  to 
meet  some  timing  criterion — for  example,  each  tank  company  arrives  at 
a specified  position  along  its  route  at  a specified  time. 

The  def ending  fire  units  need  not  be  all  of  the  same  type.  They 
might  consist  of  a mixture  of  tanks,  Sheridans,  Tows,  and  Dragons 
positioned  at  locations  in  an  Integrated  defensive  posture.  Different 
tactical  procedures  may  be  represented  in  the  model.  For  example,  each 
defense  unit  may  remain  in  position  or  it  may  fire  from  its  position 
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and  then  move  to  a nearby  position  to  avoid  return  fire.  Each  weapon 
might  have  a specific  number  of  rounds  available.  After  the  time  taken 
to  complete  the  move,  the  unit  resumes  fire  as  targets  become  visible. 

If  two  or  more  targets  become  visible  simultaneously,  the  target  with 
the  smallest  number  of  previously  assigned  rounds  is  chosen.  This  pro- 
cedure approximates  an  integrated  fire-coordination  capability  and,  in 
effect,  also  serves  as  a proxy  for  removing  targets  that  would  have 
been  destroyed  by  previous  rounds. 

The  process  of  firing,  moving,  and  occupying  a new  defense  position 
is  continued  until  the  rounds  available  are  exhausted,  until  some  tank 
company  approaches  so  close  to  the  defense  position  th.it  the  defense 
position  is  abandoned,  or  until  an  enemy  tank  company  reaches  its  ob- 
jective. 

Table  2 presents  the  computer  input  to  the  ser -Icing  rate  program 
for  this  illustration.  The  maximum  and  minimum  ranges  given  in  the 
table  are  not  the  real  range  limits  for  existing  weapon  systems.  The 
maximum  ranges  were  selected  to  represent  bad  weather  visibility  con- 
ditions and  the  minimum  ranges  were  selected  to  represent  the  range  at 
which  the  defense  units  move  out  of  their  positions  to  avoid  becoming 
overrun . 

The  output  from  the  program,  which  is  listed  in  App.  B,  is  a 
minute-by-minute  indication  of  events.  Table  3 presents  a sample  oul- 
put  at  Minute  8 of  a computer  run  involving  companies  ol  tanks  attacking 


Table  2 

EXAMPLE  OF  INPUTS  TO  SERVICING  RATE  ROUTINE 


Velocity  » 18  km/hr 

Maximum  rounds  per  target  - 4 


Input 


Weapon  1 


Weapon  2 


Weapon  3 Weapon  4 


No.  of  rounds  available  40 

Min  firing  range  750  m 

Max  firing  range  1000  ra 

Firing  rate  1 per  30  sec 

Firing  interval  60  sec 

Movement  interval  120  sec 


29 

750  m 
1000  m 
45  sec 
90  sec 
120  sec 


10 

750  m 
1000  m 
45  sec 
90  sec 
120  sec 


6 

750  m 
1000  m 
45  sec 
90  sec 
120  sec 


| 

l 

l| 

\ 


f Ai 
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Table  3 

EXAMPLE  SERVICING  RATE  ROUTINE  OUTPUT 


Cumulative  results:  Time  = 8 Minutes 


Route 

Attacker  Tank 

Number 

No. 

Name 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Total 

1 

Cl 

4 

2 

2 

2 

2 

0 

0 

0 

0 

0 

12 

2 

C2 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

3 

3 

C3 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

2 

4 

C4 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

5 

C5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

C6 

2 

2 

0 

1 

1 

0 

0 

0 

0 

0 

6 

7 

C7 

4 

1 

1 

1 

1 

l 

0 

0 

0 

0 

9 

8 

C8 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

*) 

0 

CO 

l 

1 

2 

0 

0 

0 

0 

0 

0 

0 

4 

10 

CIO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Number  ol  targets  serviced  through  this  minute  = 27 

Position  Routes 

No.  Name  Cl  C2  C3  C4  C5  C6  C7  C8  CO  CIO  Tot.  Cum. 

11  F71S  0000010000  11 

20  F72S  0000020  0 00  2 3 

21  F73S  00000100  0 0 14 

22  F74S  0000020000  2 b 

23  F75S  0000002000  2 8 

24  F7bS  0 0 0 0 0 0 1 1 0 0 2 10 

25  F81S  000  00000  1 0 111 

38  K44T  0 2 0 0 0 0 0 0 0 0 2 13 

30  E43T  2 0 0 0 0 0 0 0 0 0 2 13 

60  EblT  2 0 0 0 0 0 0 0 0 0 2 17 

61  F.62T  2000000000  2 10 

62  F.63T  2000000000  2 21 

63  F.64T  2 0 0 0 0 0 0 0 0 0 2 23 

64  K63T  2 0 0 0 0 0 0 0 0 0 2 23 

65  F.01T  0120000000  3 28 

71  F22T  0002000000  2 30 

75  F31T  0 0 0 0 0 0 0 0 2 0 2 32 

7b  F32T  0000001010  2 34 

7 7 F33T  0000002000  2 36 

78  F34T  0000001  10  0 2 38 

70  F35T  0000002000  2 40 

Total  12  322069240  40 


L. 
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on  each  of  ten  routes.  Output  data  In  the  top  half  of  the  table  show 
the  number  of  rounds  that  fired  against  each  of  the  ten  tanks  on  each 
of  the  routes  through  the  eighth  minute.  On  route  C2 , for  example 
(second  row  of  the  table),  enemy  tanks  numbered  1,  4,  and  9 were  fired 
upon  once  during  the  first  eight  minutes.  The  "Total"  column  tabulates 
the  rounds  fired  on  each  route.  During  the  eight  minutes,  27  tanks 
(the  number  of  cells  in  the  body  of  the  table  with  entries  other  than 
0)  have  been  "serviced." 

The  lower  portion  of  the  table  presents  the  cumulative  number  of 
rounds  fired  from  each  defensive  position,  indicated  by  the  identifi- 
cation code  at  the  left,  at  all  the  tanks  on  each  route.  For  example, 
defense  weapon  number  71,  Code  F22T,  has  fired  two  rounds  against  tanks 
on  Route  C4.  A cross-comparison  of  the  two  parts  of  the  table  shows 
that  the  two  rounds  fired  by  this  defense  weapon  were  directed  against 
tanks  1 and  4 of  the  tank  company  on  that  route.  As  indicated,  27  tanks 
have  been  "serviced"  through  the  eighth  minute  at  least  once,  and  a 
total  of  40  rounds  have  been  dired  by  the  defense  units.  From  the  out- 
put tables  for  each  minute,  the  entire  "servicing"  picture  can  be  con- 
structed to  determine  which  weapons  fired  how  many  rounds  at  which  enemy 
vehicles  at  what  time,  taking  into  account  the  intervisibilities  of  the 
specific  terrain. 

The  detailed  data  in  Table  3 can  be  summarized  to  show  the  cumu- 
lative number  of  enemy  vehicles  serviced  over  a period  of  time  as  well 
as  the  number  serviced  each  minute.  Figure  5 presents  an  example  of 
such  a summary  for  a force  of  100  enemy  tanks  attacking  on  ten  different 
routes.  The  figure  shows  that  for  the  specific  terrain  used  in  the 
analysis,  the  intervisibility  conditions,  combined  with  the  dispositions 
and  weapon  capabilities  of  the  defense  force,  provide  this  force  with 
the  potential  to  service  all  of  the  enemy  tanks  and  that  they  are  all 
serviced  within  12  minutes  after  the  start  of  the  attack. 


Number  of  enemy  tanks  lerviced 
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V.  USES  AND  LIMITATIONS  OK  TIMER 


\< 


USES  OF  TIMER 

TIMER  was  developed  for,  and  has  been  used  in,  several  ARPA 
studies  of  advanced  svstem  concepts.  In  these  studies,  the  model 
provided  data  on  acquisition  ranges,  time  and  distance  of  target  ex- 
posure, firing  opportunities,  and  other  factors  for  direct  fire. 
TIMER  has  also  been  used  in  somt  terrain  analyses  of  forward  defense 
positions  in  U.S.  V Corps.  Data  from  TIMER  have  also  contributed  to 
several  studies  of  airborne  platforms. 

In  view  of  the  emphasis  being  placed  on  new  concepts  of  ground 
combat  that  take  advantage  of  terrain  ("terrain  tuning") , rapid  re- 
sponse ("management  of  seconds"),  and  precision  guided  munitions 
(PCMs) , TIMER  can  be  a useful  tool  for  detailed  assessment  of  the 
Interplay  among  terrain,  intervisibility  of  weapons,  and  movement 
o f un its. 


L IMITATIONS  OK  11 HER 

The  TIMER  model  uses  a digitalized  terrain  data  base  and  a series 
of  computer  programs.  The  digitalized  terrain  data  base  is  generally 
of  good  quality,  but  several  characteristics  should  be  recognized. 
Recent  clearing  of  forests  and  new  construction  in  urban  areas  have 
not  been  incorporated  in  the  data  base.  No  attempt  has  been  made  to 
determine  the  specific  magnitude  of  these  changes,  but  they  appear  to 
be  unimportant  thus  far. 

A second  characteristic  of  the  data  base  is  that  some  local  fea- 
tures— such  as  a small  clump  of  trees,  a single  building,  a small  un- 
dulation in  the  terrain,  a tree-lined  road— are  not  included  in  the 
data  base.  These  could  influence  intervisibility  in  some  specific 
cases;  however,  for  evaluations  involving  multiple  defense  positions 
and  multiple  attack  routes  or  vehicles,  their  effects  are  considerably 
reduced.  Another  characteristic  is  that  urban  areas  are  treated  as 
bounded  areas:  "Openings"  such  as  central  squares,  plazas,  wide  roads, 
and  soccer  fields  are  not  distinguished  in  the  data  base. 


' ! i 
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These  characteristics  influenced  the  degree  of  resolution  that 
could  he  achieved  in  some  of  the  specific  computations  made  in  the 
ARP  A studies,  but  the  overall  quality  of  the  data  base  makes  it  a 
unique  and  valuable  tool  for  most  intervisibility  studies. 

The  current  computer  programs  of  TIMER  have  several  limitations. 
They  do  not  include  a capability  to  determine  the  dynamic  effects  on 
reduced  intervisibility  due  to  local  weather  conditions,  smoke,  dust, 
etc.,  on  target  servicing  rates.  They  also  do  not  include  the  capa- 
bility to  vary  the  velocity  of  moving  vehicles  in  response  to  undula- 
tions of  the  terrain  or  the  surface  conditions  of  roads  and  cross- 
country areas. 

Finally,  TIMER  is  not  a two-sided  engagement  model.  It  does  not 
incorporate  hit  probabilities,  weapons  effects,  or  exchanges  of  fire. 
In  the  evaluation  method  described  in  this  series  of  reports,  such 
factors  were  handled  by  constructing  probabilistic  Monte  Carlo  models 
and  sampling  by  random  numbers  to  produce  realistic  outcomes  of  combat 
engagements. 
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Appendix  A 


THE  TIMER  PROGRAM  AND  DATA  BASE 


The  TIMER  program  takes  inputs  characterizing  an  observer  and  the 
locations  defining  the  routes  (the  routes  are  linear  between  these  lo- 
cations), and  determines  the  route  segments  (equal-length  portions  of 
the  route  starting  from  the  first  location)  that  are  visible  to  the 
observer  (see  Table  1).  Several  types  of  output  are  constructed  from 
these  basic  intervisibility  data: 

o A frequency  table  of  the  number  of  contiguous  visible 
route  segments  as  a function  of  sensor  range? 
o A frequency  table  of  the  number  of  firing  opportunities 
as  a function  of  the  velocity  of  the  target? 
o A stylized  map  with  a resolution  of  100  m and  containing 
a letter  or  other  symbol  indicating  the  location  of  the 
observer  and  a set  of  symbols  indicating  route  segments 
in  forest  or  town  (*),  the  symbol  of  the  observer  for 
visible  segments  (the  letter  or  other  symbol),  and  a 
decimal  point  if  the  segment  is  not  visible  and  not 
in  a forest  or  town; 

o An  output  "intervisibility"  data  set  containing,  for 
each  observer /route  combination,  a table  whose  place 
positions  represent  the  consecutive  segments  along  the 
route  and  whose  entries  are  the  ranges  from  the  observer 
to  the  route  segments,  if  the  segments  are  visible. 

THE  BASIC  LINE-OF-SIGHT  CALCULATION 

The  determination  of  a line  of  sight  between  two  points  on  a map 
is  conceptually  straightforward:  Simply  draw  a straight  line  between 
the  two  points  and  determine  whether  any  point  on  the  earth's  surface 
(plus  forest  height  or  urban  height  if  appropriate)  in  the  plane  deter- 
mined by  the  line  and  the  earth's  center  lies  above  the  line.  In 
TIMER,  the  following  calculational  procedure  was  implemented: 
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Let  PI  - (xl,yl,hl)  and  P2  - (x2,y2,h2)  be  the  Universal  Transverse 
Mercator  (UTM)  coordinates  and  elevation  above  the  earth’s  surface  (m) 
of  the  observer  and  target  locations,  respectively.  Let 

d - sqrt  | (xl  - x2)2  + (yl  + y2)2) 

be  the  (approximate)  distance  between  the  two  points  at  the  earth's 
surface.  Next,  let  HI  and  H2  be  the  elevation  (m)  of  the  two  locations, 
as  measured  relative  to  the  tangent  plane  at  the  observer  location,  so 
that  111  * hi  and 

H2  - h2  - 500  d 2 / K , 

where  K is  the  earth’s  radius  (km)  and  the  expression  subtracted  from 
l>2  is  the  approximate  drop  below  the  tangent  plane  due  to  the  earth's 
curvature.  Next,  divide  the  distance  d by  the  grid  size  (12.5  meters) 
to  obtain  the  number  of  equally  spaced  intermediate  points  N between 
PI  and  P2  at  which  the  line  of  sight  is  to  be  tested.  Each  intermediate 
point  is  then  tested  by  determining  its  elevation  above  the  tangent 
plane  (equal  to  its  elevation  obtained  by  interpolating  in  the  data 
base  minus  the  drop  due  to  the  earth's  curvature  between  the  observer 
location  and  the  test  point  plus  the  forest  or  urban  height  if  any) 
relative  to  the  value  obtained  by  interpolating  between  HI  and  H2  on 
distance  between  PI  and  the  test  point.  If  the  tost  elevation  is  larger 
than  the  interpolated  value,  the  line  of  sight  is  blocked  at  that  point; 
if  it  is  smaller,  test  the  next  point. 

ORGANIZATION  OF  THE  TIMER  DATA  BASE 

The  accurate  determination  of  a line  of  sight  between  two  points 
on  or  near  the  surface  of  the  earth  requires  a fine-grained  representa- 
tion of  the  topography  intervening  between  the  two  points.  The  Topogra- 
phic Center  of  the  Defense  Mapping  Agency  has  developed  such  a fine- 
grained data  base  for  selected  portions  of  Western  Europe,  primarily 
in  the  Federal  Republic  of  Germany.  The  DMA  data  base  consists  of  the 
elevation  in  meters  and  an  indicator  variable  (=■  1 if  the  point  is  in 
forest;  * 2 if  in  a town;  » 0 otherwise)  for  each  data  point  on  a 12.5 
m grid  for  a number  of  rectangular  regions  in  latitude  and  longitude; 
each  rectangle  encompasses  an  area  12'  in  latitude  by  24 ' in  longitude — 
i.e.,  50°  12'  to  50°  24'  in  latitude  and  9°  12’  to  9°  3b’  longitude. 
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The  coordinate  origin  for  each  rectangular  region  is  the  Intersection 
of  the  southernmost  latitude  and  westernmost  longitude,  and  the  12.5  m 
grid  is  in  12.5  m increments  in  UTM  coordinates.  Northing  and  Easting, 
starting  from  the  origin. 

To  obtain  a data  base  for  an  area  contained  in  the  DMA  data  base 
regions,  the  portions  needed  must  be  combined  and  t lie  various  origins 
kept  track  of,  or  data  points  must  be  interpolated  for  the  grid  points 
desired.  For  use  in  applications  of  TIMER,  data  bases  for  several 
sub-areas  included  within  the  overall  DMA  data  base  were  prepared.  The 
data  base  most  often  used,  and  referred  to  as  the  TIMER  data  base,  is 
a portion  in  the  border  area  between  West  and  East  Germany  in  the  U.S. 

V'  Corps  sector.  A rectangular  area  ( in  UTO  coordinates)  approximately 
2000  km" is  represented,  the  region  included  between  5592-5654  km 
Northing  and  547-579  km  Easting  in  UTM  zone  12. 

The  large  number  of  data  points  (the  data  value  at  each  data  point 
is  equal  to  ten  times  the  elevation  plus  the  value  0,  1,  or  2 of  the 
indicator  variable)  in  the  TIMER  data  base  (3+  million)  and  the  large 
number  of  elevation  evaluations  required  to  determine  the  existence 
of  a line  of  sight  between  an  observation  point  and  a set  of  points 
defining  a path  or  route,  necessitated  the  adoption  of  a computa- 
tionally efficient  packaging  ol  the  data  base  for  retrieval  of  eleva- 
tion values.  The  TIMER  data  base  grid  of  12.5  m was  arranged  In  1 km 
by  1 km  blocks  of  1600  data  points,  starting  with  the  southwest  corner 
of  the  data  and  proceeding  to  the  northwest  corner  and  repeating  from 
west  to  east. 

The  resulting  data  base  of  62  (5b54-5592  km  Northing)  x 32  (579- 
547  km  Easting)  = 1984  blocks  of  1600  data  points  (or  3,174,400  total 
data  points)  were  put  on  an  IBM  2314  disk  pack  for  permanent  storage. 

The  TIMER  program  was  designed  so  that  the  data  base  permanent  storage 
disk  pack  could  be  used  directly,  or  any  rectangular  subregion  of  the 
TIMER  data  base  could  be  written  onto  the  IBM  system  disk  packs  and  used 
by  the  program.  (The  latter  procedure  shortens  turnaround  time  because 
the  IBM  2314  disk  pack  does  not  have  to  be  mounted  for  each  computer 
run.)  The  computer  program  was  written  so  that  three  of  the  1 km  by 
1 km  blocks  of  data  points  were  always  in  the  computer  core  memory: 
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the  block  containing  the  observation  point,  the  block  containing  the 
current  route  point,  and  an  intermediate  block  along  the  line  of  sight 
containing  the  current  point  at  which  elevation  is  to  be  determined 
(if  not  in  the  other  two  blocks).  This  procedure  is  a cost-effective 
tradeoff  between  computer  core  for  data  point  storage  and; input  opera- 
tions for  reading  in  the  data  points. 

V 

TIMER  COMPUTATIONAL  SEQUENCE 

The  TIMER  program  first  reads  an  input  card  defining  the,  map  cen- 
ter and  the  number  of  5 km  by  12  km  maps,  in  a north-south  map  strip, 
to  be  printed  showing  the  routes.  Then  the  program  reads  the  data  set 
defining  the  routes.  The  last  "route"  is  not  actually  a route  ^ut  is 
used  to  indicate  the  path  of  the  border  between  the  two  Germanies  on 
the  "maps"  that  are  printed  containing  the  routes;  this  last  route  may 
be  used  to  represent  any  other  continuous  descriptive  feature.  The 
maps  of  the  routes  are  then  printed. 

Next,  three  cards  are  read  defining  parameter  values  that  will  be 
common  to  all  the  observers:  the  route  segment  length,  the  minimum 
visibility  range,  the  number  of  map  sections  in  the  map  indicating  the 
observer/route  intervisibility,  the  total  number  of  observers,  11  maxi- 
mum visibility  ranges  for  defining  the  table  of  intervisibility  versus 
range,  and  eight  velocity  values  for  defining  the  table  of  visibility 
time  intervals  versus  velocity. 

A card  is  then  read  for  each  observer  containing  the  observer's 
coordinates,  a four-character  name  for  identifying  the  observer  on  out- 
put listings  and  maps,  the  heights  of  the  observer  and  target,  the 
heights  of  the  forests  and  towns,  and  two  angles  defining  the  observer's 
coverage  sector. 

Starting  with  the  first  segment  of  the  first  route,  the  program 
determines  the  existence  or  nonexistence  of  a line  of  sight  to  each 
consecutive  route  segment.  As  the  process  proceeds,  the  program  con- 
structs an  output  table  of  first  and  last  members  of  stretches  of 
visible  segments  as  a function  of  range  (see  Table  1).  Using  this 
table  as  a basis,  two  additional  tables  are  constructed  when  the  pro- 
gram reaches  the  end  of  a route:  a frequency  table  of  the  number  of 
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stretches  of  visible  segments  as  a function  of  range,  and  a frequency 
table  of  intervisible  time  intervals  as  a function  of  velocity  and 
range.  The  latter  table  can  be  interpreted  as  the  number  of 
firing  opportunities  given  the  firing  rate.  After  all  the  routes  are 
processed,  the  program  creates  an  output  "intervisibility"  data  set 
containing,  for  each  observer/route  combination,  a one-dimensional 
table  whose  place  positions  represent  the  consecutive  segments  along 
the  route  and  whose  entries  are  the  ranges  from  the  observer  to  the 
segments  if  the  segments  are  visible,  and  zero  otherwise.  The  program 
prints  a stylized  map  with  a resolution  of  100  meters  and  containing 
a letter  or  other  symbol  indicating  the  location  of  the  observer  and  a 
set  of  symbols  Indicating  route  segments  in  forest  or  town(*),  the 
observer  symbol  for  visible  route  segments,  and  a decimal  point  if 
the  segment  is  not  visible  and  not  in  a forest  or  town. 

TIMER  INPUT  FORMATS  AND  PROGRAM  LISTING 


Tables  A.l,  A. 2,  and  A. 3 present  the  TIMER  input  formats,  the 
input  format  for  the  route  data,  and  the  program  listing,  respectively. 
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Card 

Columns 

Card  1 

1-5 

6-10 

11-15 

16-20 

21-25 

26-30 

31-35 

36-40 

Card  2 

1-10 

11-20 

21- 30 
31-40 

Card  3 

1-5 

51-55 
Card  4 
1-5 

36-40 

Card  5 

1-10 

11-20 

22- 25 
26-30 
31-35 
35-40 
41-45 
46-50 
51-55 
56-60 

Card  6 

etc. 


Table  A.l 

INPUT  FORMATS— CARD  DATA 


Data  Definitions 


Number  of  map  sections  to  be  used  for  the  map  of  the  routes 
X-coordinate  (Easting)  of  the  map  center  - 500  (km) 
Y-coordinate  (Northing)  ol'  the  map  center  - 5500  (km) 

Minimum  X-coordinate  (Easting)  of  the  data  base  - 500  (km) 
Maximum  X-coordinate  (Easting)  of  the  data  base  - 500  (km) 
Minimum  Y-coordinate  (Northing)  of  the  data  base  - 5500  (km) 
Maximum  Y-coordinate  (Northing)  of  the  data  base  - 5500  (km) 
Crid  size  (km) 


Route  segment  length  (m) 

Minimum  visibility  rnage  (m) 

Number  of  map  sections  for  map  for  each  observer 
Number  of  observers  (for  sizing  output  data  set) 


Maximum  visibility  range  (km)  number  1 
Maximum  visibility  range  (km)  number  11 

Velocity  (km/hr)  number  1 
Velocity  (km/hr)  number  8 


X-coordinate  (Easting)  - 500  (km)  for  the  observer  location 
Y-coordinate  (Northing)  - 5500  (km)  for  the  observer  location 
Four-character  name  for  the  observer 
Height  of  the  observer  above  the  surface  (m) 

Viewing  angle  number  1 (degrees  clockwise  from  north) 

Viewing  angle  number  2 (degrees  clockwise  from  north) 

Height  of  the  target  (m) 

Height  of  the  forests  (m) 

Height  of  the  towns  (m) 

Number  of  the  observer  (for  position  in  output  data  set) 


Same  as  Card  5 for  additional  observers 


Tab  le  A . 3 


TIMER  FORTRAN  LISTING 


C MAIN 

COMMON/ A/DPATH , RANGE (11) ,NSEG(25 ,2) , XROUTK ( 2000 ) , Y ROUTE ( 2000 ) , 

1 DROUTE(2000)  ,NQUAD(25 , 2 ) , IQUADX(5000) , 1QUADY  (5000) , 

2 NORS , NVQ 1 , N VQ , VQUADX  ( 1 000 ) , VQUADY  (1000),  VTQUA1)  ( 1 000 ) , 

3 1 ROUTE ( 100) , DLONG ( 100)  ,T1ME(  100)  ,VE1.(  100) , 

4 HOBS , HTARG , HFOR , HTOWN , NROUTE , NMAP , R2M l N , R2MAX , 

5 SYMBOL  ( 25 , 4 ) , NAME  ( 4 ) , XLOW , XH 1 OH , YLOW , YH 1 Gil , GS 1 ZE , NYPMAX 
1 NTF.GF.R*2  1 QUADX , 1 QllADY , JQIIADX , JQIJADY , TQUAD , 1 ROUTE 

1 NTEGER* 2 VQUADX , VQUADY , VTQUAD , NAME , SYMBOL 
1 NTEGER* 2 LENGTH ( l 6) , FTIME ( 1 6 ) 

1 NTEGER* 2 MARKR(1600) 

DATA  FT  IMF./ 5 ,10,20,30, 40 , 60 ,75,100,1 20 ,150,1 80 , 240 , 300 , 360 , 480 , 

* 600/ 

C 

C TIMER  PROGRAM  FOR  A PART  OF  GERMANY 

C NORTHING:  592-654  KM  ; EASTING:  47-79  KM  : GS1ZE=0.0125 

C 

INTEGER* 2 IF 
DATA  TF/'*7 

DEFINE  FILE  3(  2 70, 3200,1... IK) 

C 

NOBS=0 

NVQ=0 

C 

C READ  IN  ROUTES  AND  LIST  (LAST  'ROUTE'  IS  ACTUALLY  BORDER) 

C 

READ (5,98)  NMAP , XCENT , YCENT , XLOW , XII 1 GH  , Y LOW , YH  1 GH , GS  I ZF. 

98  FORMAT ( 1 5 , 7 F5 . 0 ) 

IF(XLOW.EQ.O.O)  XLOW  = 47.0 
IF(XHIGH.EQ.O.O)  XH 1 GH  = 79.0 
IF(YLOW.EQ.O.O)  YLOW  = 92.0 
lF(YHlGH.EQ.O.O)  Y1I I GH  = 154.0 
IF(GSIZE.RQ.O.O)  GS1ZE  = 0.0125 

WR 1 TE  (t> , 99  ) NMAP , XCENT , YCENT , XLOW , XII I GH , Y LOW , YH  I GH  , GS  1 ZE 

99  FORMAT (' 1' , 'NUMBER  OF  MAP  PAGES  FOR  ROUTES:  ’ .13/ 

* 'O', 'MAP  CENTER  COORDINATES:  X = ’,F(>.2,';  Y = ',F6.2/ 

* 'O', 'INPUT  DATA  REGION:  EASTING  = ',2F6.1,'  NORTHING  = ’,2F(>.1, 

* ' GS I ZE  = ' , F6 . 4 , ' KM’) 

CALL  ROUTES 

CALL  MAP (XCENT , YCENT ,0,0) 

C NR  IS  THE  NUMBER  OF  ROUTES 
NR  = NROUTE  - 1 
C 

C READ  CASE  PARAMETERS 


non 
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C 

READ(5 , 100)  DPATH , RMIN , NMAP , NVPMAX , 

* RANGE, (VEL(I) ,1=1,8) 

100  FORMAT (2F1 0.0, 21 10/ 11F5.0/8F5.0) 

WRITE(6, 101)  DPATH, RMIN, NMAP, NVPMAX, RANGE, (VEL(I) ,1=1,8) 

101  FORMAT ( ' 1 ' , 20X , ' CASE  PARAMETERS '//' 0 ROUTE  STEP  SIZE  =' , 

* F5.1,'  METERS'/ ’ O' ,' MINIMUM  DETECTION  RANGE:  ',F10.0,'  METERS'/ 

* 'O', 'NUMBER  OF  MAP  PAGES  FOR  EACH  OBSERVATION  POINT: ',14/ 

* ’O’, 'TOTAL  OBSERVATION  POINTS  = ’,13/ 

* ’O’ , 'MAXIMUM  DETECTION  RANGES: ' , 1 1 ( 1X.F6.2) , ’ (KM)’/’0’, 

* 'VELOCITIES  FOR  EXPOSURE  DISTRIBUTION: ' ,8(1X,F5. 1) , ' (KM/HR)') 

C 

R2MIN=RMIN*RMIN/ (1000*1000) 

READ  IN  OBSERVER  LOCATIONS 

1 READ(5 , 102 ,END=50)  PX, PY, NAME, HOBS, ANGLE 1 ,ANGLE2,HTARG,HFOR, 

* HTOWN.NVP 

102  FORMAT (2F 10 . 0 , 1X,4A1 ,6F5 .0,110) 

NOBS=NOBS+l 

C DETERMINE  ELEVATION  OF  OBSERVATION  POINT 
CALL  ELE V ( PX , PY , 1 , ALT , TOPOG ) 

C PLACE  OBSERVATION  POINT  SYMBOL  IN  MAP  TABLE 
NVQ  = NVQ  +1 
NVQ1  = NVQ 

I =10.0*(PX-XLOW+0.05)+1 .001 
J =10.0*(PY-YLOW+0.05)+1 .001 
VQUADX(NVQ)  = I 
VQUADY(NVQ)  = J 
VTQUAD(NVQ)  = NAME(4) 

R2MAX=RANGE (11 )**2 

C DETERMINE  OBSERVER  SECTOR  OUTLINE  FOR  MAP 
IF ( ANGLE2-ANGLE 1 . NE . 360 . 0 ) 

* CALL  QUADS(PX,PY,ANGLE1 ,ANGLE2,NVQ) 

WR I TE ( 6 , 1 04 ) PX , PY , ALT , TOPOG , NAME , HOBS , ANGLE 1 , ANGLE2 , HTARG , 

* HFOR,HTOWN 

104  FORMATC 1' .'OBSERVATION  POINT:  XCOORD  =',F6.2,'  YCOORD  =',F6.2, 

* ' ALTITUDE  =' ,F7 . 0 , ' TOPOG  =' ,F3 . 0 ,’ SYMBOL  = ’ ,4A1/ 

* ' ',' OBSERVER . HT . =',F5.0, 

* ' V. ANGLE  1 =' ,F6.0, ' V. ANGLE  2 =’,F6.0, 

* ' TARG.HT  =' ,F6. 1 , ' FOR.HT  =' ,F5. 1 , ' TOWN  HT  = ’ ,F5 . 1 ) 
WRITE(6 , 105)  (FTIME(I),I=1,16) 

105  FORMATCO’ ,54X, 'REQUIRED  PATH  LENGTH  (METERS)’// 

1 ' ' ,52X, ' FIRING  TIME  (SEC.)'/ 

2 ' ' ,25X, ' VEL. ' ,714,915/'  ' ,25X, ' (KM/HR) ' ) 

DO  3 1=1,8 

DO  2 J=1 , 16 

2 LENGTH (J)  = (VEL(I)*FTIME(J))/3.6  + 0.5 

3 WRITE(6, 106)  VEL(I ) , (LENGTH(J) , J=1 , 16) 

106  FORMAT ( ' ' , 25X.F5 . 0 , IX, 714 ,915) 
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C 

C DETERMINE  INTERV1S1B1L1TY  AND  NUMBER  OF  FIRING  OPPORTUNITIES 
C FOR  CURRENT  OBSERVER 
C 

DO  10  N = 1 ,NR 

C ZERO  OUT  OUTPUT  INTERVISIBILITY  ARRAY 
DO  4 J=1 ,lb00 

4 MARKR  (,J)=0 

NS  1 = NSEG(N,1) 

NS2=  NSEG(N,2) 

C TEST  FOR  LINE  OF  SIGHT  TO  SEGMENT  NS  OF  ROUTE  N 
DO  5 NS=NS1 ,NS2 

R2  = IPX-XR0UTE(NS))**2  + (.PY-YROUTE(NS)  )**2 
C TEST  IF  SEGMENT  IS  WITHIN  RANGE  AND  WITHIN  VIEWING  ANGLE 
IFIR2.GT.R2MAX)  GO  TO  5 

CALL  V ANGLE INOBS , PX , PY  , ANGLE  1 , ANGLES , XROUTE (. NS ) , YROUTElNSl , INOUT) 
1F( I NOUT . EQ. 1 ) GO  TO  b 

5 CONTINUE 
GO  TO  7 

b WRITE  (.6, 107)  PX , PY  , NAME , ( SYMBOL  I N , K ) , K- 1 , 4 ) 

1 0 7 FORMAT ( ' 1' , ' OBSERVAT ION  POINT:  X = ' , Fb . 2 , ' Y = ' , Fb . 2 , ’ SYMBOL= ’ , 

* 4A1 , ' ROUTE  : ' ,4A1) 

C DETERMINE  FIRING  OPPORTUNITIES  ON  ROUTE  N FOR  CURRENT  OBSERVER 
CALL  FIRE (N , PX , PY , ANGLE  1 , ANGLE 2 . MARKR ) 

7 1 F ( NVP . F.Q . 0 . OR . NVPMAX . EQ . 0 ) GO  TO  10 
JK=NVPMAX*(N-1)+NVP 
WRITE(3’JK)  MARKR 
10  CONTINUE 

C PRINT  MAP  SHOWING  VISIBLE  SEGMENTS  ON  ALL  ROUTES 

WRlTE(b, 104)  PX, PY, ALT, TOPOG, NAME, HOBS, ANGLE  1 .ANGLES .HTARG, 

* HFOR , HTOWN 

CALL  MAP ( PX , PY , 1 , 0 ) 

GO  TO  1 
50  CALL  EXIT 
END 
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SUBROUT INK  F 1 RK (NP , PX , PY , ANGLE  1 , ANGLE2 , MARKR ) 

COMMON/ A/DPATH .RANGE ( 1 1 ) ,NSEG(25 ,2) ,XR011TE(2000) ,YR0UTE(2000) , 

1 OROUTE (2000) ,NQUAD( 25,2) , 1QUADX(5000) , 1QUADY(5000) . 

2 NOBS  ,NVQl , NVQ , VQUADX  (.1000),  VQUADY  (1000),  VTQUAD  ( 1 000 ) , 

.1  1 ROUTE 1 100)  ,DLONG(  100)  ,TIME(  100)  ,VF.L(  100) , 

4 HOBS , HTARG , HFOR , HTOWN , N ROUTE , NMAP , R2M 1 N , R2MAX , 

5 SYMBOL  (25,4),  NAME  ( 4 ) , X LOW , XH I GH , Y LOW , YH  l GH , GS 1 ZF. , NVPMAX 
1 NTKGER*2  I QUADX , 1 QUADY , JQUADX , JQUADY , TQUAD , 1 ROUTE 

1 NTECER* 2 VQUADX , VQUADY , VTQUAD , NAME , SYMBOL 
1 NTEGER*2  MARK ( 1 1 ) , F 1 R 1NG ( 1 1 , 1 bO ) , I FI RE( 3 1 ) , TF , MARKR ( 1 ) 

DATA  TV'/'*’/ 

DATA  RADIAN/0.0174533/ 

C NP=PATH  NUMBER;  NS=SEGMENT  NUMBER;  N=TF.ST  POINT  NUMBER 
WRITE (6, 100)  RANGE 

1 00  FORMAT ( ' 0 ' , 6 1 X , ' F 1 R l NG  OPPORTUN 1 TY  SEGMENTS ' / 

* ' ROUTE’ ,63X,' MAX  RANGE  (KM)’/ 

* ' '.’POINT  XCOORD  YCOORD  RANGE  BEARING', 

* ’ ' , 11(F8.2)/) 

DO  1 1=1,11 
MARK ( I ) = 0 

DO  10  J=l , 1O0 
1 FIRING! I, J)  = 0 
I FLAG=0 
MARKP  = 1 
DCUM  =0.0 
N POINT  = 1 
NS  = NSEG(NP.l) 

X = XROUTE(NS) 

Y = YROUTE(NS) 

5 R2  = (X-PX)**2  ♦ (Y-PY)**2 
INVIKW=1 

1 F ( R2 . GT . R2MAX ) 1NVIEW=0 
IF(R2 . LT.R2M1N)  1NVIEW=0 
IF(INVIEW.EQ.O)  GO  TO  10 

CALL  VANGLE(NVQ1,PX,PY,ANCLE1 ,ANGLE2,X,Y, INVIEW) 

IF(INVIEW.EQ.O)  GO  TO  10 

CALL  LOS  ( NPO 1 NT , PX , PY  , X , Y , XT , YT , HTOPOG , 1 NV I F.W ) 

I F ( 1 NV I EW . EQ . 0 . AND . (XT . NE . X . OR . YT . NE . Y ) ) GO  TO  10 
C SET  MAP  SYMBOL 

I = 10*(X-XLOW+0.05)  ♦ 1.00 l 
J = l 0* ( Y -YLOW+O . 05 ) ♦ 1.001 

II  = VQUADX ( NVQ ) 

JJ  = VQUADY (NVQ) 

lF(l.EQ.II.AND.J.EQ.JJ)  GO  TO  10 
NVQ  = NVQ  ♦ 1 

IF(NVQ.GT.  1000)  WRITK(6,U7) 

117  FORMAT ('O' , ' THE  NUMBER  OF  VIEWED  QUADS  EXCEEDS  THE  MAX  OF  1000') 
1F(NVQ.GT  1000)  CALL  EXIT 
VQUADX (NVQ)  = I 
VQUADY (NVQ)  = J 


VTQUAD(NVQ)  = NAME(4) 

1F( INVIEW. EQ. 1)  GO  TO  10 
VTQUAD(NVQ)  = TF 
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C 

C PRINT  FIRING  OPPORTUNITY  SEGMENTS  AND  SET  FIRING  OPPORTUNITY  TABLE 
C 

10  IF(1NVIEW.EQ.1. AND. NP01NT.LE. 1600)  MARKR(NPOINT)= 1000*SQRT(R2) 

I SEE  = 12 

IF(INVIEW.EQ.O)  GO  TO  12 
DO  11  1=1,11 

11  I F( RANGE (12-1 )**2 . GT. R2 ) ISEE=12-I 

C 1SF.E  IS  NO.  OF  SMALLEST  RANGE  LARGER  THAN  DISTANCE  TO  TEST  POINT 
IF(ISEE.EQ.l)  GO  TO  13 
IF (MARK ( I SEE- 1 ) . EQ. 0)  GO  TO  13 

12  I F l MARK (ll).EQ.O)  GO  TO  13 

C PRINT  PREVIOUS  TEST  POINT  AS  F.ND  OF  SEGMENT  IF  NOT  ALREADY  PRINTED 
IF(MARKP.EQ. 1)  GO  TO  13 
R = SQRT(R2L) 

C R2L  IS  SQUARED  RANGE  TO  PREVIOUS  TEST  POINT 
TEMPI  = XL  - PX 
TEMP2  = YL  - PY 

BEAR  = 90.0  - AT AN2(TEMP2, TEMPI) /RAD IAN 
NN=NPOINT- 1 

WR  l TF.  ( 6 , 1 0 1 ) NN ,XL,YL,R,BEAR,  (MARK(I)  ,1  = 1,11) 

C SET  FIRING  OPPORTUNITY  TABLE 

13  DO  16  1=1,11 

I F ( I . LT . I SEE)  GO  TO  14 
MARK!  I)  = MARK(  0 + 1 
C TEST  FOR  LAST  TEST  POINT  ON  ROUTE 

l F ( NS+ 1 . EQ . NSEG (NP , 2 ) . AND . DCUM+DPATH/ 1 000 . 0 . GT . DROUTE (NS+ 1 ) ) 

* GO  TO  14 
GO  TO  16 

14  M = MARK(I) 

IF(M.F.Q.O)  GO  TO  16 
M=M1N0(M , 160) 

FIRING(I.M)  = FIRING (I, M)  + 1 

IFLAG=1 

MARK ( I ) = 0 

16  CONTINUE 

C MARK  CURRENT  TEST  POINT  AS  NOT  PRINTED;  IE  NOT  INITIAL  SEGMENT  POINT 
MARKP  = 0 

IF(  I SEE . F.Q.  12)  GO  TO  20 
I F (.MARK  ( I SEE  ) . NE  . 1 ) GO  TO  20 

C PRINT  CURRENT  TEST  POINT  AS  INITIAL  POINT  OF  A SEGMENT 

17  R = SQRT(R2) 

TEMPI  = X - PX 
TEMP2  = Y - PY 

BEAR  = 90.0  - ATAN2(TEMP2,TEMP1)/RADIAN 
WRITE (6,101)  NPO I NT , X , Y , R , BEAR , ( MARK (I), 1=1, 11) 

101  FORMAT ( ’ ' ,I5,2F7.2,F6.1,F7.1,1X,UI8) 
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C MARK  CURRKNT  TEST  POINT  AS  PRINTED 
MARKP  = l 

20  DCIIM  = DCUM  ♦ DPATH/ 1000.0 

23  1F(DCUM. LE ,DR0UTE(NS+1 ))  GO  TO  24 
NS  = NS+1 

IF(NS.GE.NSEG(NP,2))  GO  TO  30 
GO  TO  23 

24  NP01NT  = N POINT* 1 
XL=X 

Y1.=Y 

R2L=R2 

DEL  =(DCUM-UROUTElNS) )/ (DROUTElNS* 1 )-l)ROUTE(NS) ) 

X = XROUTEINS)  + DE  L*  (,  XROUTE  ( NS  + 1 ) - XROllTE  (NS)) 

Y = YROUTE(NS)  + DEL*(YROUTE(NS+  l)-YROUTE(NS) ) 

GO  TO  S 
C 

C PRINT  FIRING  OPPORTUNITY  TABLE 

C 

30  1 F ( 1 FLAG . EQ . 0 ) RETURN 
WRITE (6, 102)  RANGE 

102  FORMAT L ’ 0 ' , 29X , ' PATH ' , 27X , ' NUMBER  OF  FIRING  OPPORTUNITIES'/ 

* ’ ’ , 28X . ’ LENGTH ' , 35X , ' MAX  RANGE  ( KM ) ' / 
l ’ ’,2bX.  ’ METERS V ,11(F8. 2) /) 

DO  40  J=1 , lbO 
JJ  = DPATH* J 
DO  4 l 1=1,11 
IFIREl 1 1=0 
1X1  JK=.J,  1 bO 

IS  IFIREl l 1=1FIRE( l 1+FIR1NG( 1 , JKl*(JK/J) 

I F 1 1 F I RE  1 1 1 1 . F.Q . 0 1 GO  TO  SO 
WRITE  lb, 1031  JJ, l IFIREl 1 1 .1  = 1,111 

1 03  FORMAT l ’ ’ , 2SX . 1 8 , 1 1 1 8 1 
40  CONTINUE 

C 

C PRINT  EXPOSURE  11  ME  DISTRIBUTION 


SO  WRITE  lb,  1 101  PX , PY  .NAME  , l SYMBOL  l NP , 1 1 , 1 = 1 .-^l 

110  FORMAT  1’ 1 ',  'OBSERVATION  POINT:  X =',Fb.2,’  Y =’,Fb.2,'  SYMBOl.=  ’ , 

* 4A 1 , ' ROUTE  : ' , 4A 1 / ' 0 ’ , 4SX , ' EXPOSURE  T l ME  D 1 S I R 1 BU  I- 1 ON  ' 

* ' ' ,4SX, ' (20  SECOND  INTERVALS!' /, 


•>v 

iV 


'0','  VEL  RANGE  0 20  40  bO  80  10*12 
’ 28*30* 32*34* 3b*38*40*42*44*4b*48*S0*S 
’ TOTAL  AVE  MAX'/ 


14* lb*  1 8*20*22*  24*  2 b* ' . 
*S 4*Sb*S 8* OVE R lOlAl  ' . 


’ KM/HR  KM  20  40  bO  80  10* 12* 14* lb* 18*20*22*. 
’ 30* 32*34* 3b* 38* 40* 42*44* 4b* 48*S0*S2*S4*Sb* S8*b0*b00 


* ’ ( SEC)  (SEC)  (SEC)’) 

DO  bO  IVEL= 1 ,8 

IF(VEL(  IVF.L)  .EQ.O.O)  GO  TO  bO 

1 F 1 1 VEL . EQ . S ) WR 1 TE  l b , 1 1 0 ) PX . l'Y  , N AME , ( SYMBOL  (NP , l ) , 1 = 1 ,4) 
WRlTE(b, 1 11) 


8*  ’ 


i 

i 

i 

f 


i 

I 

■j 

i 
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111  FORMAT ( ' ' ) 

DO  60  IR=1 ,11 
XTIMET=0 . 0 
NTOTAL=0 
MAXT=0 

DO  54  K=1 ,31 

54  IFIRE(K)=0 

DO  55  J=1,160 

XTIME=J*DPATH*3.6/VEL(IVEL) 

K=XTIME/20.0+l 
IF(K.GT. 30)  K=31 
NFIRE=FIRING(IR, J) 

IFIRE (K)=IFIRE(K)+NFIRE 
XT  I MET=XT  I MET+XT I ME*NF  I RE 
NTOTAL=NTOTAL+NF I RE 
IF(NFIRE.NE.O)  MAXT=XTIME+0 . 5 
ITIMET=XTIMET+0 . 5 

55  CONTINUE 

KAVE=XT I MET/ MAXO ( 1 , NTOTAL ) +0 . 5 

WRITE(6, 112)  VEL(IVEL) ,RANGE(IR) , (IFIRE(K) ,K=1 ,31) .NTOTAL, 
* ITIMET , KAVE ,MAXT 

112  FORMAT ( 1 ’ ,F5. 1 ,F6. 2, 2X, 3013, 14,416) 

60  CONTINUE 

RETURN 

END 
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SUBROUTINE  LOS (POINT, XI ,Y1 ,X2,Y2 ,XF ,YF,HTOPOG,LOS) 

COMMON/ A/DPATH, RANGE ( 11) ,NSEG(25,2) ,XR0UTE(2000) ,YROUTE(2000) , 

1 DROUTE (2000), NQUAD ( 25 , 2 ) , I QUADX ( 5000 ) , I QUADY (5000), 

2 N0BS,NVQ1 ,NVQ,VQUADX(1000) ,VQUADY(1000) ,VTQUAD(1000) , 

3 IROUTE(IOO) ,DLONG( 100) , TIME (100) ,VEL(100) , 

4 HOBS , HTARG , HFOR , HTOWN , NROUTE , NMAP , R2MIN , R2MAX , 

5 SYMBOL (25,4), NAME ( 4 ) , XLOW , XHI GH , YLOW , YH I GH , GS I ZE , NVPMAX 
INTEGER*2  I QUADX , I QUADY , JQUADX , JQUADY ,TQUAD , IROUTE 
INTEGER*2  VQUADX , VQUADY , VTQUAD , NAME , SYMBOL 

REAL*8  X,Y,H,DELX,DELY,DELH 

C LOS  = 1 INDICATES  LINE  OF  SIGHT  BETWEEN  X1,Y1  AND  X2,Y2 
C LOS  = 0 INDICATES  NO  LINE  OF  SIGHT 

C XF,YF  IS  THE  FIRST  TEST  POINT  AT  WHICH  LINE  OF  SIGHT  IS  OBSCURED 
C HTOPOG  IS  TOPOG  AT  XF,YF 
DATA  RE/6378.388/ 

LOS  = 1 

DIST  = SQRT((X1-X2)**2+(Y1-Y2)**2) 

IF(DIST.LE.GSIZE)  RETURN 
CALL  ELEV (X2 , Y2 , 2 ,H2 , TOPOG) 

H2  = H2  + HTARG 

IF (TOPOG . EQ . 0 . 0)  GO  TO  1 

IF (TOPOG. EQ. 1.0. AND. HTARG. GT. HFOR)  GO  TO  1 

IF (TOPOG. EQ. 2.0. AND. HTARG. GT. HTOWN)  GO  TO  1 

LOS  = 0 

X = X2 

Y = Y2 
D = 0.0 
H = 0.0 
HACT  = 0.0 
HTOPOG  = HFOR 

IF (TOPOG. EQ. 2.0)  HTOPOG  = HTOWN 
GO  TO  6 

1 CALL  ELEV (XI ,Y1,1, HI, TOPOG) 

HI  = HI  + HOBS 
N = DIST/GSIZE 
DELX  = (X2-X1)/(N+1) 

DELY  = (Y2-Y1)/ (N+l) 

C CORRECT  H2  FOR  EARTH  CURVATURE 

DELH  = (H2- 1000. 0*0. 5*DIST*DIST/RE  - HI)/ (N+l) 

DELD  = DIST/ (N+l) 

X = XI 

Y = Y1 
H = HI 
D = 0.0 

C 

DO  5 M=1,N 
X = X + DELX 

Y = Y + DELY 
H = H + DELH 
D = D + DELD 
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CALL  ELEV(X,Y,3,HACT,T0P0G) 

C CORRECT  HACT  FOR  F.ARTH  CURVATURE  AND  TEST  FOR  LINE  OF  SIGHT 
HACT  = HACT  - 1000. 0*0. 5*D*D/RE 
ITOPOG  = TOPOG 
HTOPOG=0 . 0 

I F( ITOPOG. EQ.O)  GO  TO  4 
GO  TO  (2,3,2),  ITOPOG 

2 HTOPOG  = HFOR 
GO  TO  4 

3 HTOPOG  = HTOWN 

4 I F ( H . GT . HACT+HTOPOG ) GO  TO  S 
LOS  = 0 

GO  TO  6 

5 CONTINUE 
C 

6 XF=X 
YF=Y 

IF(HACT.GT.H)  HTOPOG=O.Q 
999  RETURN 
END 
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SUBROUT  1 NK  ELEV  (X , Y , 1 N , ALT , TOPOG ) 

COMMON/ A/DPATH .RANGE ( 1 1 ) ,NSEG(25,2) , XROUTE ( 2000 ) .YROUTE(2000) , 

1 DROUTE(2000) ,NQl)AD(25,2) , IQUADX(5000) , IQUADY(5000) , 

2 NOBS , NVQ 1 , NVQ , VQUADX  (1000),  VQIIADY  (1000),  VTQ11A1)  (1000), 

3 1 ROUTE (100), DEONG ( 1 00 ) , TI ME ( 1 00 ) , VEL ( 1 00 ) , 

4 HOBS , HTARG , HFOR , HTOWN , NROUTE , NMAP , R2M 1 N , R2MAX , 

5 SYMBOL (25,4),  NAME  ( 4 ) , XLOW , XH 1 Gtl , YLOW , YH I GH , GS 1 ZE , NVPMAX 
1 NTEGER* 2 1 QUADX , l QUADY , JQUADX , JQUADY , TQllAl) , 1 ROUTE 

1 NTEGER* 2 VQUADX, VQUADY .VTQUAD , NAME , SYMBOL 
C 

C ELEVATION  IN  METERS  AND  TOPOGRAPHY  FOR  POINTS  IN  A SELECTED  PART  OF 
C GERMANY.  TOPOGRAPHY:  1=F0REST;  2=URRAN ; 0-NE1THER  FOREST  NOR  URBAN 
C NORTHING:  592-654  KM  ; EASTING:  47-79  KM 

C X = EASTING (KM) ; Y = NORTHING (KM) ; ALT  = ELEVATION (M) 

C 


1  NTEGER'V2  Z (40 , 40 , 3) , Z 1 ( 1 600 ) , Z2 ( 1600 ) , Z3(  1600) 

EQUIVALENCE  (Z(  1 , 1 , 1 ) ,Z1 ( 1 ) ) , (Z( 1 , 1 ,2) , Z2( 1 ) ) , (Z( 1 , 1 ,3) ,Z3( 1 ) ) 
INTEGER  H ( 2 , 2 ) , I H ( 2 ) , 1JG(3) 

REAL  HY ( 2 ) , XORG ( 8 ) , YORG ( 8 3 
DATA  XORG/  0.002,  0.009,  0.004,  0.010, 

fc  0.011,  0.011,  0.010,  0.008/ 

DATA  YORG/  0.004,  0.006,  0.008,  0.012, 

* 0.010,  0.011,  0.001,  0.004/ 


DATA  lJG/3*0/ 

DEFINE  FILE  1 (7936, 3200, L, 1J) 


C 

C THE  PROGRAM  FROM  HERE  TO  FORTRAN  STATEMENT  12  IS  AN  AD  HOC  WAY  TO 
C FIND  THE  ORIGIN  CORRECTIONS  FOR  THE  VARIOUS  PORTIONS  OF  THE  DATA  BASE 
C IT  WOULD  NOT  BE  NEEDED  FOR  OTHER  DATA  BASES 


IF(Y.LT.YL0W+0. 0125. OR. Y.GT.YH1CH-0. 0125)  GO  TO  13 
lF(X.LT.XL0W+0. 0125. OR. X.CT.XHIGH-O. 0125)  GO  TO  13 
IF(X.GE. 71.0)  GO  TO  1 
IF(X. LE. 70.0)  GO  TO  2 

D = 0.0135254*(Y-628. 176)  + 0 . 9999085*(X-70 . 473) 

I F( ABS (D) . LE . 0 . 0 1 25 ) GO  TO  12 
IF(D.LT.O.O)  GO  TO  2 

1 D = (0 . 999868* (Y-628 . 1 76) +0 . 0 1 62 28* (X- 70 . 473) )/22 . 242  + 7.0 
GO  TO  3 

2 D = (0 . 999942* (Y-628 . 1 76)+0 . 0 1 0799*(X-70 .473))/22. 242  + 3.0 

3 I D=D 

IF(D-ID.LE. 0.001 124)  GO  TO  12 
IF(D-ID.GE. 0.998876)  GO  TO  12 
XOR=  XLOW+XORG( I D) 

YOR=  YLOW+YORG( ID) 

IFLAG=0 
GO  TO  20 
12  1 Fl,AG=  1 

XOR=  XLOW+GS I ZE 
YOR=  YLOW+GS I ZE 
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GO  TO  20 
13  ALT=0 . 0 
TOPOG=0 . 0 
RETURN 

20  XI  = (X-XOR)/GSIZE 

I = XI 

YJ  = (Y-YOR)/GSIZE 
J = YJ 
1)X  = XI  - I 
DY  = YJ-J 
DO  30  11=1,2 
DO  29  JJ=1 ,2 
1ND  = IN 

IG  = (I+II-D/40 
JG  = (J+JJ-D/40 
IJ  = 2*(YHIGH-YL0W)*IG+JG+1 
DO  21  K=1 ,3 

IF(IJ.NE. IJG(K))  GO  TO  21 
IND  = K 
GO  TO  26 

21  CONTINUE 

GO  TO  (22,23,24),  IND 

22  READ(l'lJ)  Z1 
GO  TO  25 

23  READ(l'IJ)  Z2 
GO  TO  25 

24  READ(l'IJ)  Z3 

25  lJG(IND)  = IJ-1 

26  H(ll ,JJ)  = Z(l+II-40*IG,J+JJ-40*JG,IND) 
I F ( I FLAG . EQ . 1 ) GO  TO  40 

29  IH(JJ)  = H(II ,JJ)/10 

HY(II)  = (IH(2)-IH(1))*DY  + IH(1) 

30  CONTINUE 

ALT  = (HY(2)-HY(1))*DX  + HY(1) 

I = 2*DX  + 1 
J = 2*DY  + 1 

TOPOG  = H(I , J)  - 10*(H(I , J)/ 10) 

RETURN 

40  ALT  = H( 1 , 1 ) / 10 

TOPOG  = H ( 1 , 1 ) - 1 0*ALT 

RETURN 

END 
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SUBROUT I NE  MAP (XCENT , YCENT , MODE , MT I ME ) 

COMMON/ A/ DPATH , RANGE ( 1 1 ) , NSEG ( 25 , 2 ) , XROUTE ( 2000 ) , YROUTE ( 2000 ) , 

1 DROUTE(2000) ,NQUAD(25 ,2) , IQUADX(5000) , IQUADY(SOOO) , 

2 NOBS , NVQ 1 , NVQ , VQUADX (1000), VQUADY (1000), VTQUAD (1000), 

3 I ROUTE ( 100) ,DLONG(100) ,TIME(100) ,VEL(100) , 

4 HOBS , HTARG , HFOR , HTOWN , NROUTE , NMAP , R2M 1 N , R2MAX , 

5 SYMBOL (25,4), NAME ( 4 ) , XLOW , XH I GH , YLOW , YH I GH , GS I ZE , NVPMAX 
I NTEGER*2  I QUADX , I QUADY , JQUADX , JQUADY , TQUAD , I ROUTE 
INTEGER*2  VQUADX , VQUADY , VTQUAD , NAME , SYMBOL 

DIMENSION  HEAD (10) 

INTEGER" 2 A( 12 1 ) , B (26 ) , P ,TF ,OH , BORDER , BLANK 
DATA  TF/'*7 

DATA  P/ * . OH/'O’/, BORDER/' Z ’ / , BLANK/'  '/ 

DATA  HEAD/0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9/ 

C 

C PRINTS  MAPS  CENTERED  AT  XCENT, YCENT 
C 

IMIN  = 1 0*A I NT ( XCENT - 5 . 0 - XLOW ) + 1 
JMIN  = 1 0 " A I NT ( YCENT - ( 5*NMAP ) / 2 - YLOW ) + 1 

11  = XCENT- 5 . 0 

12  = 11+12 
DO  1 1=11,12 
11=1-11+1 
B(II)=I/10 
B(II+13)=I-10*(I/10) 

1 CONTINUE 
C 

IF(MTIME.NE.O)  WRITE(6,100)  MTIME 
IF(MTIME.EQ.O. AND. MODE. EQ.O)  WRITE(6,104) 

100  FORMATC  1' , 'TIME  =',I4) 

WRITE(6 , 105)  B 
NLINE=50*NMAP+ 1 
DO  20  J=1 ,NLINE 
JJ  = NLINE-J+1 
DO  11  1=1,121 

11  A ( I ) = BLANK 
DO  12  N=l, NROUTE 
K1  =NQUAD(N , 1 ) 

K2  =NQUAD(N,2) 

DO  12  K=K  1 , K2 

IF ( IQUADY (K) . NE . JJ+JMIN- 1 ) GO  TO  12 
IF(IQUADX(K).LT.IMIN)  GO  TO  12 
IF(IQUADX(K) .GT. IMIN+120)  GO  TO  12 
II  = IQUADX(K)-IMIN+1 
A(II)  = P 

IF (MTIME. EQ.O. AND. MODE. EQ.O)  A(II)  = SYMB0L(N,4) 

IF(N.EQ. NROUTE)  A(II)  = BORDER 

12  CONTINUE 

14  IF (NVQ. EQ.O)  GO  TO  17 


...  


C 
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M1  = NVQ1 
M2  = NVQ 

IF(MODE.EQ.O)  Ml=l 
DO  16  M=M1  ,M2 

IF(VQUADY (M) . NE . JJ+JMIN- 1 ) GO  TO  16 
IF(VQUADX(M)  . LT.  1MIN)  GO  TO  16 
IF(VQUADX(M).GT. IMIN+120)  GO  TO  16 
II=VQDADXIM)-IMIN+1 
IF(MTIME.EQ.O)  GO  TO  15 
IF(AUl).EQ.P)  GO  TO  16 
1F(VTQUAD(M) .EQ.TF)  GO  TO  16 
A ( 1 1 ) = OH 
GO  TO  16 

15  A ( 1 1 ) = VTQ'JADIM) 

16  CONTINUE 

17  JJJ  = MOD(JJ-l.lO)  +■  1 

JH  = AlNT(YCENT-(5*NMAP)/2)  + (JJ-l)/lO 
YH=JH 

H = HEAD (JJJ) 

IF(JJJ.EQ.l)  WKITE(6 , 101 ) Yll.A.YH 

101  FORMAT ( ' ' ,F6. 1 , 1 2 1 A 1 ,F5. 1 ) 

IF(JJJ.NE.l)  WRITE (6 , 102)  H , A , H 

102  FORMAT ( ' ' ,4X, F2 . 1 , 121 Al , 3X, F2 . 1 ) 

IF(JJ.EQ.l)  WRITE(6 , 103)  R 

103  FORMAT  ('  ’,6X,12('  \2345b78<) ' )/ ' ' ,bX,  12(11,' '),U/ 

* ’ ’ ,6X, 12(1 1 ,9X) ,11) 

I F(MOD( JJ ,50) . NE . 2 . OR . JJ . EQ. 2)  GO  TO  20 
WRITE (6, 103)  B 
WRITE (6, 104) 

104  FORMAT! ’17) 

WRITE (6, 105)  B 

1 05  FORMAT ( ’ ’ ,6X,  12(11 ,9X)  ,11/,’  \6X,  12(11,’ * ) . 1 1/ 

* ’ ’ , 6X , 1 2 ( ’ 123456789’)) 

20  CONTINUE 

RETURN 

END 
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SUBK0UT1NE  QUADS (XI ,Y1 ,X2 ,Y2 ,CUMDIS , J) 

COMMON/ A/DPATH , RANGE ( 1 1 ) , NSKG(2*3 ,2) ,XROUTE(2000) ,YR0UTE(2000) , 

1 DROUTE ( 2000 ) ,NQUAD(25,2) , 1QUADX(5000) , 1QUA1)Y(5000) , 

2 NOBS , NVQ l , NVQ , VQUADX (1000), VQUADY (1000), VTQUAD (1000), 

3 1 ROUTE  (100),  DLONG ( 100 ) , T 1ME  ( 1 00)  , VE1. ( 1 00 ) , 

4 HOBS , IITARO , HKOR , HTOWN , N ROUTE , NMAP , R2M 1 N , R2MAX , 

5 SYMBOL (25 , 4 ) , NAME ( 4 ) , XLOW, XH 1 GH , YLOW , YH 1 GH , GS 1 ZE , NVPMAX 
INTEGER* 2 1QUADX, IQUADY , JQUADX , JQUADY , TQUAD , IROUTE 

1 NTEGER*2  VQUADX , VQUADY , VTQUAD , NAME , SYMBOL 
C DETERMINES  QUADS  ENTERED  BY  THE  LINE  SEGMENT  FROM  XI, Y1  TO  X2,Y2 
C FOR  USE  WITH  MAP 

C CUMDIS  IS  THE  CUMULATIVE  DISTANCE  ALONG  THE  ROUTE 
C J IS  CURRENT  NUMBER  FOR  QUAD  TABLE 
X =10. 0*(Xl-X1.0W+0 . 05 )♦  1 . 00 1 
Y =10. 0* ( Y l -YLOW+O . 05 ) ♦ 1 .001 
DIST  = 10 . 0*CUMDIS 
IX  = X 
1Y  = Y 

IF(CUMDIS.NE.O.O)  GO  TO  10 
J=J+1 

IF(J.GT.5000)  WRITE (b, 100) 

1F(J . GT .5000)  CALL  EXIT 
1QUADX(J)  = IX 
IQUADY (J)  = IY 
10  DELX  = tO*(X2-Xl) 

DF.LY  = 10*(Y2-Y  1 ) 

D = SQRT  (DELX**  2 + DF.LY**  2) 

IF(DEl.X.EQ.O.O)  TX  = 100.0 
1 F (DELX. EQ. 0.0)  GO  TO  1 
INCX  = 1 

DE1.TX  = ABS(  1 . 0/DEl.X) 

FRACX  = X - IX 

TX  = (1.0  - FRACX) /DELX 

IF(DELX.GT.O.O)  GO  TO  1 

INCX  = -l 

TX  = -FRACX/ DELX 

1 IF(DELY.EQ.O.O)  TY  = 100.0 
IF(DELY.EQ.O.O)  GO  TO  2 
INCY  = 1 

DELTY  = ARS( t .O/DELY) 

FRACY  = Y - IY 

TY  = (1.0  - FRACY ) /DF.I.Y 

IF(DELY.GT.O.O)  GO  TO  2 

INCY  = -1 

TY  = - FRACY / DEL Y 

2 T = AMIN1 (TX,TY) 

IF(T.GT. 1.0)  GO  TO  5 
J = J + 1 

1F(J  .LE.5000)  GO  TO  b 
WR1TE(6,100) 


100  FORMAT ('O' , ' THE  NUMBER  OF  QUADS  EXCEEDS  THE  MAX  OF  5000') 
CALL  EXIT 

6 IF(TX.GT.TY  + 0.001)  GO  TO  3 
IX  = IX  + INCX 
TX  = TX  +DELTX 

IF(TY.LT.TX  -DELTX  + 0.001)  GO  TO  3 
GO  TO  4 

3 IY  = IY  + INCY 
TY  = TY  + DELTY 

4 IQUADX(J)  = IX 
IQUADY(J)  = IY 
GO  TO  2 

5 CUMDIS  = CUMDIS  + D/10.0 
RETURN 


n o n a 
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SUBROUTINE  ROUTES 

COMMON/ A/DPATH , RANGE (11), NSEG (25,2), XROUTE ( 2000 ) , YROUTE ( 2000 ) , 

1 DR0UTE(2000) ,NQUAD(25 ,2) , IQUADX(5000) , 1QUADY(5000) , 

2 NOBS , N VQ 1 , NVQ , VQUADX ( 1 000 ) , VQU ADY ( 1 000 ) , VTQUAD (1000), 

3 I ROUTE ( 100) ,DLONG( 100) ,TIME(100) ,VEL(100) , 

4 HOBS , HTARG , RFOR , HTOWN , NROUTE , NMAP , R2M I N , R2MAX , 

5 SYMBOL(25 ,4) .NAME (4) ,XLOW,XHIGH,YLOW, YHIGH,GSIZE,NVPMAX 
I NTEGER* 2 1 QUADX , I QUADY , JQUADX , JQUADY , TQUAD , I ROUTE 
INTEGER*2  VQUADX , VQU ADY .VTQUAD .NAME , SYMBOL 
I NTEGER* 2 SYM(4) ,STEMP(4) 

READS  IN  ROUTE  SEGMENTS  AND  LISTS 

ALSO  DETERMINES  QUADS  ENTERED  BY  THE  ROUTES  (FOR  MAP) 

NROUTE  = 0 
DO  20  1=1,4 
20  STEMP( I )=0 . 0 
C READ  IN  ROUTE  SEGMENTS 
J=0 

NSEG( 1 , 1)=1 

1 READ(2, 100,END=10)  SYM,X,Y 

100  FORMAT (6X,4At ,2Ft0.2) 

J = J + 1 

IF(J.GT.2000)  WR1TE(6, 101) 

101  FORMAT ('O' , ' THE  NUMBER  OF  ROUTE  POINTS  EXCEEDS  THE  MAX  OF  2000') 
IF(J.GT. 2000)  CALL  EXIT 

XROUTE (J)  = X 
YROUTE(J)  = Y 
DO  2 1=1,4 

IF(SYM( 1 ) . NE . STEMP( I ))  GO  TO  3 

2 CONTINUE 

NSEG ( NROUTE , 2 ) =NSEG ( NROUTE , 2 ) + 1 
GO  TO  l 

3 NROUTE  = NROUTE + 1 
DO  4 1=1,4 

SYMBOL (NROUTE, I)  = SYM(I) 

4 STEMPll)  = SYM(I) 

l F ( NROUTE . NE . 1 ) NSEG ( NROUTE , 1 ) =NSEG ( NROUTE - 1 , 2 ) + 1 
NSEG ( NROUTE , 2 ) =NSEG ( NROUTE , 1 ) 

GO  TO  1 

C PRINT  ROUTE  SEGMENTS 

10  WRITE(6 , 102) 

102  FORMATC 1’ ,55X, 'ROUTES’/ ' O' , ’ROUTE' ,5X,5('  SEG.  XCOORD  YCOORD ' ) ) 
DO  12  N=l, NROUTE 

WRITE(6, 103) 

NN  = NSEG(N , 1) 

NNN  = MIN0(NSEG(N,1)+4,NSEG(N,2)) 

11  WRITE(6 , 103)  (SYMBOL(N.I) ,1=1 ,4) , 

* (M, XROUTE (M), YROUTE (M),M=NN, NNN) 

103  FORMAT ( * ' ,4A1 ,6X,5(I7,F7.2,F7.2)) 
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NN  = NN+5 

NNN  = MIN0(NNN+5 ,NSEG(N ,2) ) 

I F (NN . LE . NSEG (N , 2 ) ) GO  TO  11 
12  CONTINUE 

C DETERMINE  QUADS  ENTERED  BY  ROUTES  (FOR  MAP) 
K=0 

NQUAD( 1 , 1)=1 
DO  15  N=1 ,NROUTE 

IF(N.NE. 1)  NQUAD(N,1)=NQUAD(N-1,2)+1 
NQUAD(N ,2)=NQUAD(N , 1 ) 

JJ  = NSEG(N, 1 ) 

JJJ  = NSEG(N ,2) 

DROUTE(JJ)  =0.0 
CUMD1S  =0.0 
DO  15  J=JJ,JJJ 
X = XROUTE(J) 

Y = YROUTE(J) 

IF(J.EQ.JJ)  GO  TO  15 
X1=XR0UTE(J-1) 

Y1=YR0UTE(J- 1 ) 

CALL  QUADS(X1,Y1,X,Y,CUMDIS,K) 

DROUTE ( J ) =CUMD I S 
NQUAD(N ,2)=K 
15  CONTINUE 
REWIND  2 
RETURN 
END 
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SUBKOUT 1 N E V ANGLE ( N , X , Y , ANGLE 1 , ANGLE 2 , XT , YT , 1 NOUT ) 

C DETERMINES  IF  POINT  XT , YT  IS  WITHIN  VIEWING  ANGLE  FROM  X,Y 
C VIEWING  ANGLE  IS  FROM  ANGLE  1 TO  ANGLE2  IN  CLOCKWISE  DIRECTION 
C BOTH  ANGLES  ARE  MEASURED  FROM  NORTH  DIRECTION  ON  GRID 
C INOUT  = 0 IF  POINT  XT , YT  IS  NOT  WITHIN  VIEWING  ANGLE 
C INOUT  = 1 IF  POINT  XT , YT  IS  WITHIN  TH  VIEWING  ANGLE 
DATA  RAD 1 AN/ 0.01 74533/ ,NP/0/ 

IF(NP.F.Q.N)  GO  TO  5 

NP=N 

PX=X 

PY=Y 

A1=ANGLE1 

A2=ANGLE2 

I F ( A 1 . GE . A2 ) A2 =A2 + 3b0 . 0 
C SET  UP  VECTORS  FOR  VIEWING  ANGLE 
XI  = COS( (90-A1 )*RADIAN) 

Y 1 = SIN( (90-A1  )'VRADIAN) 

X2  = COS(  (90-A2)‘vRADlAN) 

Y2  = SINK  90-A2 ) •' RAD  l AN ) 

C DETERMINE  IF  POINT  IS  WITHIN  VIEWING  ANGLE  (BY  PARALLEL  PROJECTIONS) 
5 I F ( A2  - A 1 . NE . 3b<> . 0 ) GO  TO  1 
1N0UT=1 
RETURN 

1 TEMPI  = (XT-PX)-'Yl  - (YT-PY)*Xl 
TEMP2  = (XT-PX)*Y2  - (YT-PY)*X2 
TEMP)  = X1*Y2  - Y1*X2 
IF(180.0-(A2-A1))  2,3.4 

2 INOUT  = 0 

l F ( TEMP 1/TEMP3.GT.Q.0)  1 NOUT= 1 
1 F ( TF.MP2 / TF.MP3 . LT  .0.0)  1 NOUT=  1 
RETURN 

3 INOUT  = 1 

IF (TEMP1.LT. 0.0)  INOUT  = 0 
RETURN 

4 INOUT  = I 

1 F ( TEMP  1 /TEMP 3 . GT . 0 . 0 ) 1 NOUT=0 
I F(TEMP2/TF.MP3 . LT  .0.0)  1 NOUT=0 
RETURN 
END 
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Append  lx  B 

SERVICING-RATE  PROGRAM  EXAMPLE 

The  TIMER  program  of  App.  A creates  an  output  intervisibility 
file  containing,  for  each  observer/route  combination,  a table  whose 
place  positions  represent  consecutive  equal-length  segments  (usually 
25  m)  along  the  route  and  whose  entries  are  the  ranges  from  the 
observer  to  the  centers  of  the  route  segments  if  the  segments  are 
visible  from  the  observer  and  zero  otherwise. 

This  appendix  contains  a description  and  program  listing  for  a 
computer  routine  that  uses  the  TIMER  output  intervisibility  file  as  an 
input.  This  model  calculates  the  servicing  rate  (see  Sec.  IV)  for 
firing  units  that  are  in  defensive  positions  at  the  observer  locations 
and  companies  of  enemy  tanks  approaching  along  the  routes. 

Table  B.l  presents  the  input  formats  and  Table  B.2  presents  an 
example  of  the  servicing-rate  program  FORTRAN  listing. 
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Table  B.l 
INPUT  FORMATS 


Card 

Columns 

Card  1 

1-10 

11-20 

21-30 

Card  2 

1-10 

11-20 

Card  3 

1-10 

11-20 

21-30 

31-40 

Card  4 

1-10 

11-20 

21-30 

31-40 

Card  5 

Card  6 

1-10 

11-20 

21-30 

31-40 

Card  7 

1-10 

11-20 

21-30 

31-40 

Card  8 

1-10 

11-20 

21-30 

31-40 


Data  Definitions 


Number  of  observers 
Number  of  routes 
Number  of  route  segments 


Tank  company  velocity  (km/hr) 

Maximum  number  of  rounds  fired  at  each  tank 


Number  of  rounds  for  firing  unit  type  1 
Number  of  rounds  for  firing  unit  type  2 
Number  of  rounds  for  firing  unit  type  3 
Number  of  rounds  for  firing  unit  type  4 


Minimum  range  (km)  for  firing  unit  type 
Minimum  range  (km)  for  firing  unit  type 
Minimum  range  (km)  for  firing  unit  type 
Minimum  range  (km)  for  firing  unit  type 

Same  as  card  4 for  maximum  range 


Time  (sec) 
Time  (sec) 
Time  (sec) 
Time  (sec) 


to  fire  one  round 
to  fire  one  round 
to  fire  one  round 
to  fire  one  round 


for  firing 
for  firing 
for  firing 
for  firing 


Firing  interval 
Firing  interval 
Firing  interval 
Firing  interval 


(sec)  before  moving  for 
(sec)  before  moving  for 
(sec)  before  moving  for 
(sec)  before  moving  for 


Movement 

Movement 

Movement 

Movement 


interval 

interval 

interval 

interval 


(sec)  for 
(sec)  for 
(sec)  for 
(sec)  for 


firing  unit 
firing  unit 
firing  unit 
firing  unit 


1 

2 

3 

4 


unit  type  1 
unit  type  2 
unit  type  3 
unit  type  4 


ii 

firing  unit  type  1 
firing  unit  type  2 
firing  unit  type  3 
firing  unit  type  4 


type 

type 

type 

type 


1 

2 

3 

4 
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Card 

Columns  Data  Definitions 

Card  9 

1-5  Starting  segment  for  tank  company  on  route  1 

71-75  Starting  segment  for  tank  company  on  route  15 

Card  10 

1-1  Type  (1  to  4)  of  firing  unit  at  observer  location  1 

80-80  Type  (1  to  4)  of  firing  unit  at  observer  location  80 

Card  11 

1-1  Type  (1  to  4)  of  firing  unit  at  observer  location  81 

22-22  Type  (1  to  4)  of  firing  unit  at  observer  location  102 

Card  12 

1-4  Four  character  name  for  route  1 


57-60  Four  character  name  for  route  15 
Card  13 

1-4  Four  character  name  for  observer  1 


80-80  Four  character  name  for  observer  20 
Card  14-18 

Continue  names  for  observers  21  to  102 
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Table  B.2 

SERVtCING-KATK  PROGRAM  FORTRAN  LISTING  EXAMPLE 


C MAIN 

I NTEGER*  2 M I NKGK  (4)  , MAXKGK  ( 4 ) . FRATE  ( 4 ) . FT  IMF.  ( 4 ) , MOVE  ( 4 ) . ROUNDS  ( 4 ) 
INTEGER-'.'  START  ( 102) ,FIRST(  102) .LAST ( 102)  .SHOT  ( 102) 

INTEGERS  TYPE  ( 102) , I TOTAL 1 102) ,F( 102, IS) ,FF( 102 , IS) 

1 NTEC.KR* 2 l N I T 1 1 S ) . G 1 1 S . 1 0 ) . GG  ( 1 S . 1 0 ) . P ( 1 S . 1 0 ) , .ITOTAL(  1 S ) 

I NTEGER  X ( 1 02 , 1 S , 2 7 ) , MASK ( 12 ) 

INTEGER* 2 MARK (800) 

DATA  MASK/ 7 I , 72 . 7.4 . 78 .710, 720 . 740 , 780 . 7 1 00 , 7200 . 7400 . 7800 . 7 1 000 . 

1 72000 , 74000 , 78000 , 7 1 0000 . 720000 , 740000 . 780000 , 7 1 00000 , 7200000 . 

2 7400000 . 7800000 , 7 1 000000 . 7.2000000 , 74000000 . 78000000 . 7 1 0000000 , 

1 720000000 , 740000000 , 780000000/ 

DIMENSION  A(  IS).  11(102) 

DATA  A/ IS*'  '/ 

C 

G EXAMPLE  SERVICING  KATE  PROGRAM 
C 

READ(S  , 100)  NI’OS , NROUTE , NSEG 
*00  FORMAT ( 1110) 

READ(  S .101)  VF.I. , MXSHOT . ROUNDS . M I NRGF. . MAXKGK . FKATE  , FT  I ME . MOVE 

1 0 1 FORMAT ( F 1 0 . 0 . 1 1 0/ ( 4 II 0 ) ) 

WK I I F.  (P.l  02 ) VEI.  .MXSHOT . ROUNDS . M I NRGF. . MAXRGF. , FKATE . FT  I ME . MOVE 

102  FORMAT ( ’ I ’ . 'VELOCITY  ' .FS.O, ' KM/HR ' / 

* ’O' . 'MAXIMUM  ROUNDS  PER  TARGET  '.IS / 

-•  'O'. 2 IX.'  MftOAl  '.'IX.'MSSl  '.  I0X .'  TOW'  ,8X .'  DRAGON ' / 

* 'O'. 'NO.  OF  ROUNDS:  ',4(110,'  ')/ 

* 'O’. 'MIN  FIRING  RANGE:  '.4(110.'  M ' )/ 

* ' '.’MAX  FIRING  RANGE:  '.4(110,'  M ' )/ 

'•  ' ' . ' F I R I NG  RATE : I PER  ',4(110.'  SEC ' ) / 

* ' '.'FIRING  INTERVAL:  '.4(110,'  SEC')/ 

* ’ ’.'MOVEMENT  INTERVAL: ',4(110,'  SEC')) 

READ(S,  10.1)  (INIT(N).N  I, NROUTE) 

101  FORMAT ( ISIS) 

WNITE(<>.  104)  (iNIT(N).N  I. NROUTE) 

104  FORMAT ('O',' START INC  POINTS  ON  ROUTES:  ' .ISIS) 

READ(S.IOS)  (TYPF.(N)  ,N  I ,NPOS) 

I OS  FORMAT (8011) 

WR I TK ( <» , 1 0(> ) (TYI’K(N).N  I, NCOS) 

I Oft  FORMAT ( 'O' , 'WEAPON  TYPES:  '.1021  I ) 

C A(N)  IS  A FOUR  CHARACTER  NAME  FOR  ROUTE  N 
C l»(N ) IS  A FOUR  CHARACTER  NAME  FOR  POSITION  N 
RKAIMS.lon  (A(N).N  I .NROUTE) 

RF.AIMS.IO/)  (ll(N).N  I .NCOS) 

10/  FORMAT (20A4) 

C V IS  VELOCITY  IN  METERS/ SECOND 
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V= VEL* 1000/ 3600 
C INTIAUZK  COUNTER  TABLES 

C C(I,J)  IS  CUMULATIVE  ROUNDS  ASSIGNED  TO  TANK  J OF  ROUTE  l 
C CC(I,J)  IS  CUM.  ROUNDS  ASSIGNED  TO  TANK  J ROUTE  1 THRU  PREVIOUS  MINUTE 
C F ( I , J ) IS  CUM.  ROUNDS  ASSIGNED  TO  ROUTE  J FROM  POSITION  1 
C FF(I , J ) IS  CUM.  ROUNDS  ASSIGNED  TO  ROUTE  J FROM  1 THRU  PREVIOUS  MINUTE 
DO  1 1=1 .NROUTE 
DO  1 J=l,10 
C(I,J)=0 

1 CCII,J)=0 

DO  2 1=1 ,NPOS 
DO  2 J=l, NROUTE 
F(I,J)=0 

2 FF( I , J)=0 
C 

C READ  IN  OBSERVATION  POINT/ROUTE  SEGMENT  INTERVISIBILITY  ARRAYS  AND 
C USE  BINARY  INDICATORS  FOR  INDICATING  SEGMENTS  IN  VIEW  AND  IN  RANGE 

C 

DO  5 J=l, NROUTE 

DO  5 M=l ,NPOS 

RM I N=M  1 NRGF,  ( TYPE  ( M ) ) 

RMAX=MAXRGE(TYPE(M) ) 

DO  3 K=1 , 26 

3 X(M, J ,K)=0 

READ 1 1,110)  (MARK(JK) ,JK=1 ,NSEG) 

110  FORMAT ( 20 (40A2)) 

1)0  4 K=1  ,NSF.G 

IF(MARK(K) . LE.RMIN.OR.MARK(K) .GT.RMAX)  GO  TO  4 

JJ=(K- l)/30+l 

I=K-30*(JJ-1) 

X(M, J , JJ)=X(M, J, JJ)+MASK( I ) 

4 CONTINUE 

5 CONTINUE 
C 

C INITIALIZE  EVENT  TIME  INDICATORS  FOR  EACH  OBSERVATION  POINT 
C FOR  EACH  FIRING  UNIT: 

C START  IS  THE  TIME  AT  WHICH  THE  FIRING  UNIT  WILL  NEXT  BE  READY  TO  FIRE 
C FIRST  IS  THE  TIME  AT  WHICH  THE  FIRST  ROUND  WAS  FIRED  AFTER  MOVING 
C LAST  IS  THE  TIME  AT  WHICH  THF.  MOST  RECENT  ROUND  WAS  FIRED 
C SHOT  IS  THE  NUMBER  OF  ROUNDS  REMAINING 

C 

C Nl  IS  THF.  TIME  STEP  (SEC) 

C 

N 1 = 10 

DO  6 M=1 ,NPOS 
SHOT(M)=ROUNDS(TYPE(M) ) 

START(M)=Nl 
FIRST (M)=- 1000 

6 LAST(M)=- 1000 
C 
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C REPEAT  CAECUEAT10NS  FOR  EACH  TIME  STEP  UNTIE  SOME  TANK  UNIT  REACHES 
C THE  END  OF  ITS  ROUTE 
C 

DO  90  N=N 1 , 1500 ,N 1 
NN=N/bO 

DO  10  M= 1 , NPOS 
1 F(SHOT(M) . l.E . 0)  GO  TO  30 
1TYPE=TYPE(M) 

I F ( F l RST(M ) +FT 1 ME  ( I TYPE ) - FRATE  ( 1 TYPE ) . 1.1. N) 

* START  ( M ) =F l RST ( M ) ♦ FT  I ME.  ( 1 TYPE ) +MOVE  ( I TYPE ) 

I F ( START (M ) . CiT . N ) GO  TO  30 
1 F ( EAST (M)  + FRATE ( 1 TY PE ) . GT . N ) GO  TO  30 
C DETERMINE  AEE  TANK  POSITIONS  THAT  CAN  RE  FIRED  UPON  FROM  POSITION  N 
C 1FEAG  NOT  EQUAE  TO  ZERO  INDICATES  AT  I.EAST  ONE  FIRING  OPPORTUNITY 
C FOR  THE  FIRING  UNIT  THIS  TIME  PERIOD 
I FLAG=0 

DO  15  I=l,NROUTE 
1 N I TX= 1 N 1 T (.  I ) 

DO  15  J=1 , 10 

C P(I,J)  IS  THE  CUM.  NUMBER  OF  ROUNDS  FIRED  ON  TANK  .1  OF  ROUTE  1 
C DEFAUET  VAEUE  OF  10000  INDICATES  TANK  CANNOT  BE  FIRED  UPON  THIS  PERIOD 
P(l ,J)= 10000 

1 F ( C ( I , J ) . F.Q . MXSHOT ) GO  TO  15 
lV=V*N/25-(J-l)*2  + 1N1TX 
IF(IV.LE.NSEG)  GO  TO  11 
WRITE (6 ,115) 

115  FORMAT (’ I’,'  THE  PROGRAM  HAS  RUN  OUT  OF  ROUTE  SEGMENTS') 

CAEE  EXIT 

C DETERMINE  IF  TANK  .)  OF  ROUTE  1 IS  V1S1RI.E  AND  WITHIN  RANGE  FOR  ENTIRE 
C ROUND  FIRING  TIME 

11  IF(IV.EE.  INITX)  GO  TO  15 
ID=(IV-l)/30+l 
MK=X(M,I,1D) 

IF(MK.EQ.O)  GO  TO  15 
ID=lV-30*( ID- 1 ) 

MK=MK-MASK(  ID)-MASK(  1IH1  ) )'> (MK/MASK(  1D+ 1)) 

IF(MK.ET.O)  GO  TO  15 
P(1,J)=C(I,J) 

1 W=FRATE  ( l TYPE )*V/25 . 0+0 . 5- 1 . 0 
IF(IW.EE.O)  GO  TO  14 
1X1  13  K=1 , IVV 
I v=l V- 1 

IF(IV.EK.O)  GO  TO  12 
ID=( IV- 1 )/30+ 1 
MK=X(M , I ,111) 

IF(MK.EQ.O)  GO  TO  12 
ID=IV-30,V(  ID- 1 ) 

MK=MK-MASK( ID)-MASK( ID+ 1 ) )*(MK/MASK(11)+1 ) ) 

IF(MK.NE.O)  GO  TO  13 

12  P(I,J)=10000 
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GO  TO  15 

13  CONTINUE 

14  1FLAG=1 

15  CONTINUE 

C SELECT  AS  TARGET  THAT  TANK  WITH  THE  FEWEST  NUMBER  OF  ROUNDS  ASSIGNED 
IF(IFLAG.EQ.O)  GO  TO  30 
MINP=10000 
DO  20  1=1, NROUTE 
DO  20  J=l,10 

IF(P(I,J).GE.MINP)  GO  TO  20 
MINP=P(I , J) 

II=I 

JJ=J 

20  CONTINUE 

IF(MINP.EQ. 10000)  GO  TO  30 

C(II,JJ)=C(II,JJ)+1 

F(M,II)=F(M,II)+1 

LAST(M)=N 

SHOT(M)=SHOT (M) - 1 

IF(FTRST(M) . LT. START (M))  FIRST(M)=N 
30  CONTINUE 

C 

C AT  END  OF  EACH  MINUTE  PRINT  TABLES  SHOWING  NUMBER  OK  ROUNDS  ASSIGNED 

C 

C PRINT  ROUNDS  ASSIGNED  TO  EACH  TANK  THIS  MINUTE,  BY  ROUTE 
1 F ( MOD ( N , 60 ) . NE . 0 ) GO  TO  90 
NTOTAL=0 

DO  41  1=1, NROUTE 
JTOTAL( 1 )=0 
DO  41  J=1 , 10 
CC(I,J)=C(I,J)-CC(I,J) 

JTOTAL( I )=JTOTAL( 1 )+CC( 1 ,J) 

NTOTAI.=NTOTAL+CC  ( I , J ) 

41  CONTINUE 

WRITE (6, 120)  NN, (1,1=1, 10) 

120  FORMAT  ( ' 1 ' , ' TIME  =',13,'  MIN’ ,50X, 'TARGET  //’/ 

* ' 0 ' , l 5X , 1 0 I l 0 , 1 OX , ' TOTALS ' / ) 

1X1  42  1=1, NROUTE 

42  WR ITE ( 6 , 1 2 1 ) 1 ,A( I) , (CC( I , J) , J=1 , 10) , JTOTAL( l ) 

121  FORMAT ( ' ’,15,  2X.A4.4X, 101 10, 1 15) 

WRITE(6, 122)  NTOTAL 

122  FORMAT ( ' 0 ' , 1 1 5X , 1 1 5 ) 

C PRINT  TOTAL  NUMBER  OF  TARGETS  SERVICED  THIS  MINUTE 
KTANK=0 

DO  43  1=1, NROUTE 
DO  43  J=l,10 

I F ( CC ( I , J ) . NE . 0 ) KTANK=KTANK+ 1 

43  CC(I,J)=C(1 ,J) 

C PRINT  NUMBER  OF  TANKS  SERVICED  THIS  MINUTE 
WRITE(h, 123)  KTANK 
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123  FORMAT (' O’ ,20X,  'NUMBER  OF  TARGETS  SERVICED  THIS  MINUTE  =',I4) 
C PRINT  ROUNDS  ASSIGNED  TO  EACH  ROUTE  THIS  MINUTE,  BY  POSITION 

DO  45  1=1 ,NPOS 
IT0TAL(I)=0 
DO  45  J= l , NROUTE 
FFII,J)=F(I,J)-FF(I,J) 

45  I TOTAL (I ) = I TOTAL ( I ) + FF ( I , J ) 

WRITE(b , 125)  NN, (A(I ) , 1=1 , 15) , (J, J=1 .NROUTE) 

125  FORMAT ('O’ , ' TIME  =',I3,'  MIN’ , SOX, 'ROUTES' / ' ' , 15X, 15(2X,A4)/ 

* ’ ' , 15X, 1516/) 

WRITE  (.6 , 126) 

1 26  FORMAT ( ' + ' , l 1 OX , ' TOTALS ' , 6X , ' CUM . TOTAL ’ ) 

NT0TAL=0 

DO  4b  1=1 ,NPOS 
NTOTAL=NTOTAL+ 1 TOTAL  ( 1 ) 

IF(ITOTAL(I).EQ.O)  GO  TO  4b 

WR 1 TE ( b , 127)  I ,B( l ) , (FF(I , J) ,J=1 , NROUTE) 

127  FORMAT ('  ’ , 13 ,2X,A4 ,6X, 1516) 

WR1TE(6 , 128)  I TOTAL ( I ) ,NTOTAL 

128  FORMAT ( ' + ' , 105X.I 10,1 15) 

4b  CONTINUE 

WRITE  lb, 129)  (JTOTAL(J) ,J=1 .NROUTE) 

124  FORMAT ('O' , 1 5X , 1 5 1 6 ) 

IX)  47  1=1 ,NPOS 

DO  47  J=l, NROUTE  ' 

47  FF( I , J)=F(I , J) 

C CUMULATIVE  ROUNDS  ASSIGNED  TO  EACH  TANK,  BY  ROUTE 
NT0TAL=0 

DO  51  1=1, NROUTE 
JT0TAL(l)=0 
DO  51  J=l,10 

JTOTALC 1 )=JTOTAL(  I ) +C  ( I ,J) 

NTOTAL=NTOTAL+C( 1 , J) 

51  CONTINUE 

WRITE (6, 130)  NN,(I , 1=1 , 10) 

1 30  FORMAT ('O’,'  CUMULAT l VE  RESULTS ' / ' ’ , ' TIME  = ' ,13,'  MIN' ,50X, 

* 'TARGET  //'/’O'  ,15X,  10110,1  OX,  'TOTALS'/) 

DO  52  1=1, NROUTE 

52  WR I TE ( 6 , 1 3 1 ) 1 ,A( 1 ) , (C( 1 , J) , J= 1,10) , JTOTAL( 1 ) 

131  FORMAT ('  ',15,  2X,A4,4X, 101 10, 1 15) 

WRITE (6, 132)  NTOTAL 

132  FORMAT  ('O'  , U5X.I15) 

C PRINT  CUMULATIVE  NUMBER  OF  TANKS  SERVICED 
KTANK=0 
KTANK2=0 
KTANK3=0 

IX)  53  1 = 1,  NROUTE 
DO  53  J=l,10 

I F ( C ( I , J ) . EQ . 1 ) KTANK2=KTANK2+ 1 
I F ( C ( I , J ) . GE . 2 ) KTANK2=KTANK2+2 
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IF(C(I,J).EQ.1)  KTANK3=KTANK3+ 1 
I F (C ( 1 ,J) .EQ.2)  KTANK3=KTANK3+2 
IF(C(I,J).GE.3)  KTANK3=KTANK3+3 
53  IF(C(I,J).NE.O)  KTANK=KTANK+ 1 
WR I TE ( 6 , 1 33 ) KTANK , KTANK2 , KTANK3 

133  FORMAT ('O', 20X , 'NUMBER  OF  TARGETS  SERVICED  THRU  THIS  MINUTE  =' , 

* 14/'  ', 13X, 'NUMBER  OF  (PAIRED)  TARGETS  SERVICED  THRU  THIS  MINUTE 
*=’,14/'  ' ,13X, 'NUMBER  OF  (TRIPLED)  TARGETS  SERVICED  THRU  THIS  MINI! 
*TE  ='  ,14) 

C PRINT  CUMULATIVE  ROUNDS  ASSIGNED  TO  EACH  ROUTE,  BY  POSITION 
DO  55  1=1 ,NPOS 
ITOTAL(l )=0 
DO  55  J=1 .NROUTE 

55  lTOTAL( I ) = 1T0TAL( I )+F( I , J ) 

WR1TE(6 , 135)  NN,(A(I),I=1,15),(J,J=1 .NROUTE) 

135  FORMAT ( ' 0 ' , ' CUMULATIVE  RESULTS'/'  ','  TIME  =',13,’  MIN' ,50X, 

* 'ROUTES’/'  ' , I5X, 15(2X,A4)/'  ' , 15X, 1516/) 

WR1TE(6,126) 

NTOTAL=0 
DO  56  1=1 ,NPOS 
NTOTAL=NTOTAL+ I TOTAL ( I ) 

IF(ITOTAL(I).EQ.O)  GO  TO  56 

WRITE(6, 127)  I ,B(I ) , (F(I , J) , J=1 .NROUTE) 

WRITE(6, 128)  ITOTAL(I) ,NTOTAL 

56  CONTINUE 

WRITE(6, 137)  (JTOTAL(J) ,J=1 .NROUTE) 

137  FORMAT ('O' .15X.1516) 

C 

90  CONTINUE 
END 


